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Abstract 



An Explanation of my Thesis in Words of One Syllable 

With apologies to Paul Boolos. 

My goal for the year was to learn what small things do. If you have a small thing it's not a 
point, it's a blob - we knew that from the start. But does it stay a blob? No, of course not. It 
breaks up, but at some time it comes back. At times less than that time, we find small blobs, too. 
When it's not made of small blobs, the small thing is most odd. It forms lots of hills and troughs. 
I have ways in here to find out when the blob comes back, where and when the small blobs are, 
how long it takes for the blobs to break up, where the hills and troughs are, and how many of them 
there ought to be. 

The Same Explanation with Longer Words 

The dynamics of a quantum mechanical particle in a time-independent potential are found to contain 
many interesting phenomena. These are direct consequences of the (typical) existence of more than 
one time-scale governing the problem. This gives rise to full revivals of initial wavepackets, fractional 
revivals (multiple wavepackets appearing at fractions of the revival time), and the striking quantum 
carpets. A variety of analytic techniques are used to consider the interference that gives rise to these 
phenomena while skirting calculations involving cross-terms. Novel results include a new theorem on 
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the weighting coefficients that govern fractional revivals, a demonstration that ^cl, the function 
that governs the distribution and features of these fractional revivals, really does behave classically, 
a treatment of wavepacket dephasing in the infinite square well by means of the Poisson summation 
formula, and a correct analysis of the spatial distribution of intermode traces. Also, this work 
presents a coherent treatment of these phenomena, which before now did not exist. 
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Chapter 1 



Introduction 



1.1 A Hitchhiker's Guide to the Infinite Square Well 

The infinite square well is treated in almost every introductory quantum mechanics course. It is the 
essence of simplicity, a potential bereft of features, its eigenfunctions are simple sine waves. Exact 
solutions abound. We build an intuition for how quantum mechanical objects should behave, and 
we move on - after all, the infinite square well is just too simple. 



Figure |L1| was enough to convince me that I didn't really understand the square well. I had 
no idea that its probability density was so structured - it hadn't even occurred to me that it was 
periodic. So I began the year of study that culminated in this thesis, following much the same 
pattern as before: look at the square well, find an interesting phenomenon, study it, then see if 
it occurs in other types of system^. While there are applications for most of these phenomena, I 
won't pretend that my goal was to find more of them - my goal was to gain a deeper understanding 
of a fundamental part of quantum mechanics. 



focused on time-independent potential wells, in one dimension. 
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n. Lfi 

Figure 1.1: A particle in the infinite square well. Lighter areas indicate greater probability density. 
Our initial wavefunction is a Gaussian packet, centered on x = L/2, with ax = 0.003L. For the 
definition of Tr, see Chapter |2|. 
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1.1.1 Classical Oscillation and Dispersion 
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The first interesting phenomenon is that of the classical oscillation of wavepackets. Given a well- 
localized wavepacket, or even a wavefunction with several distinct peaks, those peaks will undergo 



So long as the expectation value oi p coincides with the value of p at the center of the wavepacket, 
the packet's oscillation will be classical. Of course, this is rarely exactly true and virtually never 
stays true - when the expectation value of p and the value of p at the center of the wavepacket do 
not coincide, the wavepacket will begin to disperse, and more characteristically quantum mechanical 
phenomena will begin to appear. An obvious question to ask is, "can we quantify when this dispersion 
occurs?" A related question is "does this dispersion occur simultaneously for the whole wavefunction, 
or does it depend on which point in the wavefunction we are considering?" This question will be 
considered for a zero-dimensional "quantum beats" system in Chapter ^, and in more dimensions in 



something like classical oscillation on a time scale that will be studied in Chapter |2[ This is not 
particularly surprising if we consider Ehrenfest's Theorem^, 




(1.1) 



Chapter |. 



^See [ [Griffiths, 1995| , 17]. 
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1.1.2 Revivals 




The next phenomenon that we will consider is closely related to the previous one. It turns out 
that just as there is a time scale that governs classical oscillation, there is another that governs 
a distinctly quantum mechanical feature. At this time scale, the original configuration of the 
wavefunction reappears, in some cases approximately, in others exactly. That is, if we start with a 
localized wavepacket it will disperse, interfere, and eventually reconstruct itself - what I will call a 
revival. As part of this process, it will also undergo the same sort of initial classical oscillation that it 
undergoes immediately after its formation. Stranger still are the existence of fractional revivals, the 
appearance at times less than the revival time of multiple wavepackets, each related to the original 
wavepacket. Curiously, both full and fractional revivals became topics of study only recently, even 
though their existence should be apparent from the exp {—iEt/h) form of the time evolution terms 
of the wavefunction. This subject is the focus of Chapter and plays a part in Chapters ^ and § 
as well. 

1.1.3 Quantum Carpets 




Perhaps the most striking phenomenon is the existence of easily-discerned canals and ridges in 
the probability density, collections of which have come to be called quantum carpets. While these 
follow straight lines in the square well, they various curved lines in other potentials. These lines 
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look much like, but are distinct from, classical trajectories. In Chapter |5| we will see how the two are 
related. Carpets happen to be, I think, the most difficult phenomena to study. While the revivals 
problem, and with it the problem of classical oscillation and dispersion, can be reduced to a clever 
shuffling of phases, it is rather harder to avoid dealing with the entire parameter space of spectrum, 
eigenfunctions, and particular weighting coefficients. Considering this, it is remarkable that we are 
able to make interesting statements at all. 

1.2 An Analytic Theme: Interference Without Cross- Terms 

Last Fall I asked a friend of mine, a math major, if he knew anything about simplifying finite sums. 
"No," he replied, "I think you just have to add those up." The finite sums that I was asking about 
don't differ so much from the infinite sums that also appear in this family of problems, and they're 
ugly. I was a bit intimidated. One interesting aspect of these problems is that they all pertain to 
features that we look for after "solving" the problem. That is, we start with a set of eigenfunctions, a 
spectrum, and weighting coefficients, then try to rewrite them to make the features we are interested 
in manifest. Fortunately, we do have some analytic tools to do this. 

All of the phenomena that we are interested in are essentially interference phenomena, produced 
by the multiplication of two sums. The critical step in every technique presented here is the pack- 
aging of the interference information in a form that doesn't involve cross-terms. These techniques 
include the grouping of certain parts of the sum by phase (Chapter ^ and ^) and the application 
of a bit of analysis, the Poisson summation formula (Chapters ^ and |^. Both of these techniques 
seem initially to make things worse by introducing new sums or integrals, but those new sums and 
integrals are exactly what give us manageable forms for enough of the interference to make sense of 
things. 
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1.3 Semiclassical Physics, Distinctively Quantum Behavior 

What exactly constitutes a "semiclassical limit" is not entirely clear. Most of the time, we will take 
it to mean that the distribution of weighting coefficients n is centered around some central value, 
n, with some characteristic width An, and that the hierarchy 1 ^ An <^ n holds - i.e. that the 
packet is localized in physical space. In some cases, we will relax this to the WKB assumption that 
the potential changes much more slowly than the relevant eigenfunctions oscillate. 

It appears that all of the phenomena that I have studied appear most naturally in the study of 
semiclassical wavepackets. Although the revivals results^re quite general, they are most striking 
when we begin with a well-localized initial packet. Otherwise, we just get a mess that occasionally 
looks like the mess that we started with. Though we are able to solve the infinite square well 
exactly with each of these techniques (and there An n » is enough to produce a carpet), 
our generalizations mostly rely on the WKB approximation. Of course, since we are studying 
interference phenomena it is not surprising that our most interesting results appear when we have 
a large number of states that can interfere with each other. Nothing that we have done rules out 
the appearance of similar phenomena in superpositions of low-energy states, but we have not found 
any. 

What is remarkable is that we find so many distinctly quantum mechanical phenomena in sys- 
tems that are supposed to be becoming more like classical ones. Although we start with wavepackets 
that are more localized than would be a superposition of, say, the bottom three states in a well, 
and seem more particle-like, they exhibit bizarre behavior. What is remarkable is that, in the case 

of carpets, we can use a semiclassical tool to study this. 

^In Chapter |3| we will consider the case of quantum beats, a zero-dimensional problem that is interesting even 
when a small number of states interfere with each other. 
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I was able to talk to a lot of people, before my thesis kept me really busy. One evening I was 
talking with a friend-of-a-friend who happened to be a Chemistry major. She had been through 
Quantum Chemistry, but was still curious about what things like the uncertainty principle really 
meant. I spent a few minutes waving my hands, then went back to my room and picked up one 
of my plots of the square well. It proved to be a much more effective way of demonstrating how 
quantum mechanics can be deterministic, yet still produce so much uncertainty - this pattern is the 
particle in the box. It's obvious that it doesn't really have a velocity because it doesn't even have a 
trajectory. I think that just as quantum mechanics seems haphazard and riddled with uncertainty 
until one realizes that the wavefunction is the fundamental object in the theory, pictures seem 
ill-advised only until one realizes which objects to plot. 

I was also attracted to this subject because of a talk that included pictures. There were very 
few pictures in my otherwise superb quantum mechanics textbook, and as a consequence there were 
things that I didn't understand. It didn't bother me because I didn't know that I didn't understand 
them. The reason that this field developed so recently, and not before, is that these phenomena 
are immediately apparent in pictures and hopelessly hidden in eigenfunction expansions - and, of 
course, we couldn't generate those pictures without computers. I've included a lot of pictures in 
this thesis because I've found that while they're not worth a thousand equations, they're a fine 
substitute for a paragraph of my prose. 



1.5 A Bit of Background 



Much of this thesis is based on three fine papers: the original, authoritative article on fractional re- 
vivals, [ Averbukh and Perelman, 1989 |, a strong treatment of quantum beats, [ Leichtle et al., 1996b| , 



Leichtle et al., 1996a l, and a well-conceived (if roughly executed) paper on intermode traces, [ Kaplan et al., 200C 



Kaplan et al., 1998(] . 

There are a variety of other papers available. Most studied are revival phenomena, which have 
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been treated in a variety of ways in a variety of systems. Though the paper mentioned above is re- 



markable for its clarity and generality, other interesting works include [ Aronstein and Stroud, 1997 , 
Aronstein, 2000| , Bluhm et al., 1996, Chen and Yeazell, 1998, Jie et al., 1998, Knospe and Schmidt, 1996 



Loinaz and Newman, 1999| , |Razi Naqvi et al., 200l| , [Rozmej and Arvieu, 1998[| . A variety of ap- 
proaches to quantum carpets have been proposed, though most are restricted to the infinite square 
well. Some of the more interesting of these involve an analysis in terms of Wigner functions, 
which I have not pursued. In any case, interesting carpet papers include [ Friesch et al., 2000| , 
Grossmann et al., 1997| , [Hall et al., 2001| , [MarzoH et al., 1998b| , [Marzoh et al., 1998a|] . For a some- 



what dated discussion of some experimental aspects of these phenomena, see ||Averbukh and Perelman, 1991 1. 

One of the more intriguing applications of quantum carpet techniques is in the field of Bose- 
Einstein Condensation. Although I had little time to devote to them, the growing literature suggests 
that this is a subject of current interest: |Choi et al., 2001 , Ruostekoski et al., 2001 , Wright et al., 1997 1. 

Papers have been written which study the fractal geometry of various problems, including 
[ Berry, 1996| , [Wojcik and Zyczkowski, 2001[ , [Wojcik et al., 2000|] . The latter two focus on preparing 
wavepackets that are variations on the Weierstrass function (which everywhere continuous, nowhere 
differentiable) , while the former focuses on the fractal dimension of a particle with a smooth spatial 
distribution in an n-dimensional box. The most useful insight to come from any of this, I think, is 
that fractals proper only emerge from initial states that are not proper solutions to the Schrodinger 
equation. I still suspect that there may be some treatment of carpets that exploits their obvi- 
ous self-similarity, but the difficulty of fractal proofs combined with the difficulty of working with 
almost-fractals turned me away from this approach. 

There have been attempts to connect the carpet problem to quantum chaos, [ Provost and Baranger, 199^ , 
Saif, 2000[| . There is a sort of similarity between the phenomena of quantum carpets and chaotic 



scars (representation of particular classical orbits in the probability densities of a chaotic potential), 
though little has yet come of this. There is also one paper purporting to connect all of this to 
quantum computing, [Harter, 2001 1, though I must confess that I was unable to follow it. 

The final connection, which I would have liked to spend more time exploring, is between the prob- 
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ability density of a quantum mechanical particle in one spatial dimension and and electromagnetic 
wave propagating in two spatial dimensions, in the paraxial approximation. In that field, revivals 
are instances of the Talbot effect, and fractional revivals instances of the fractional Talbot effect. 
Apparently, the square well problem, a toy problem in quantum mechanics, translates into a wave 
guide problem of greater practical importance. For the interested reader, the primary papers on this 
subject seem to be [ Berry and Klein, 1996 , Berry and Bodenschatz, 1999 , Dubra and Ferrari, 1999| , 



Lock and Andrews, 1992 1 



Chapter 2 



Quantum Revivals, Full and Fractional 



2.1 In the Beginning 

Our quantum mechanics course began with the Schrodinger equation, 

5^ h"^ d"^^ 



and promptly separated it, assuming that the potential V depended only on x: 



<l(x.t) = ■4>(x)I{t), (23) 

inyi = E, (2.3) 

ti^ (fib 

+Vi^ = Ei,. (2.4) 



2m dx"^ 

Since we could solve Equation e]3, we did, 



f{t) = e-'T\ (2.5) 



and turned our attention to the Time-Independent Schrodinger Equation. We quickly found that 
for bound particles the possible energies E were quantized, making the general solution to the 

16 
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^{x,t) = J2^nMx)e-''^'. (2.6) 

n 

This eigenfunction expansion is useful for algebraic purposes, and admittedly we often want to 
do things with the wavefunction once we've found it. Unfortunately, the eigenfunction expansion 
conceals the various interesting things that happen during the time-evolution of a wavefunction. 
We'll start our analysis of these phenomena with the simplest, quantum revivals. 



2.2 Fun Revival Facts 

Given a well-localized distribution of weighting coefficients Cn centered around some mean value n, 
we can perform a Taylor expansion^ of En around n: 

Er,^En + Ei{n - n) + ^E'^{n - nf + ^^E'^^in - nf + .... (2.7) 



We now define 



27T Ei'^ 



(2.8) 



a set of time scales. Prom these we pick out two important time scales, a classical period T^i = Ti and 

the revival time Tr = T2 . It is the existence of more than one time scale that makes time-evolution 

in time-independent potentials interesting. 

We ignore the j > 3 terms for two reasons - they make our calculations hardeif], and in the 

semiclassical case, which we are most interested in, they really are ignorable. There are two ways 
^From here onward, I will assume that = 2m — 1. 

^If we really want to, we could consider the so-called "superrevivals." These aren't so absurd as to be unobservable, 
but the big qualitative change in behavior comes from introducing the second time scale, Tr. 



CHAPTER 2. QUANTUM REVIVALS, FULL AND FRACTIONAL 1 
to see this. First, we assume that our En can be written as a finite sum of weighted powers of n, 



M 



En — ^ d„ 



n 



(2.9) 



m=N 



This won't always be the case, but I think it's fair to consider transcendental spectra to be unusual. 
From looking at we see that we can write the jth derivative as 



M 



= Cfcn' X 



k=N 



1 



n 



m 



(2.10) 



m=k—j+l 



We know that the jth characteristic time scale {T^i is the first, is the second) will be proportional 



to 1/ 



E., 



U) 



. We also know that in the semiclassical case, n » 1. If n is larger than k, it will force the 



E^ terms toward zero, which will in turn force the higher time scales toward infinity. Another way 
to see this is to realize that the limiting spectrum of any deep potential well is at most in? . You can 
think of this as a consequence of the fact that no well can have walls harder than the infinite square 



well, but for a more detailed explanation you should consult [ Nieto and Simmons, 1979 |. Finally, 
we observe that in the semiclassical limit we are typically concerned with a small region of n-space, 
within which the first two terms of the Taylor expansion should be an adequate approximation. 

To the extent that T^ <C Tr (which is true in most cases), we can identify two distinct behaviors 
in the time-evolution of a wavepacket. Writing out our approximate wavefunctionQ, 



^'(x, t) = exp {—2TTiEnt) Cnipnix) exp 



(n — n)t (n — n)'^t 

-2TTi I h 



Tr 



(2.11) 



we can see that when t ~ T^;, the I/Tr term will make no significant contribution to the time- 
evolution, and the wavefunction will oscillate classically with period T^. With time the I/Tr term 
causes the wavepacket to disperse, and interference between the various states produces a quantum 

carpet. However, as t approaches within a few T^i of Tr, the I/Tr term contribution again becomes 

^Note that we are writing our wavefunction as a function of x rather than as a function of x. This is because 
when considering revival problems, we aren't particularly concerned about the dimensionality of the wavefunction. 
We will be concerned about this in subsequent chapters, and when we are the arrow will come off of the x. 
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small and the wavepacket returns to classical oscillation. When t = Tr, a quantum revival occurs. 
If Tr is an integer multiple of Td, the wavefunction returns to its original configuration, and if not 
then it returns to something very close to its original configuration. 



This is illustrated in Figure |2.l| , where we see the initial evolution of a particle in an infinite 



square well of width o. We use a Gaussian distribution of weighting coefficients with n = 40 and 



cn+k = ^e-^'/2-^. (2.12) 



Referring to Equation 2^ and Section |A.1.1| , we can calculate the time scales for the square well, 



^2 

Tci = — , (2.13) 

nvr 

27-2 

Tr = , (2.14) 

TT 

and find that = SOTd. We can see clearly the initial classical oscillation of the particle (the 
oscillations are straight lines, as we would expect in the square well) and its subsequent dispersion 
and interference. Note that the periods of the oscillations, so long as they are well-defined, are T^. 

Finally, we may note that in cases where the spectrum is strictly linear in n, such as the simple 
harmonic oscillator, we would (from this analysis) expect no dispersion, little to no interference, no 
quantum carpet - just classical oscillation. Though this is not entirely true, it is roughly true in 
virtually all cases^. An example of classical oscillation is shown in Figure 2^. There, n = 5 and 
an = 2 . 

1 1 1 1 2 

We also find revival-esque phenomena at t < Tr. If we examine ||^'|| at a few particular times 
in the infinite square well, as in Figure (insert figure ISW-Cuts)|2j, we notice something interesting. 
At t = Tr/2, ||^'(x,0)||^ has been refiected about the center of the well, and at t = Tr/4 there 
are two copies of the original probability distribution, one refiected about the center of the well. 

Although this particular result does not hold generally^, it may suggest to us that interesting things 

*In Chapter |5| we will learn the black art of "quadratizing spectra," which will let us make wavefunctions with 
carpets in problems with linear spectra. 

®For this result to hold we require eigenstates of definite parity and a purely quadratic spectrum. For a discussion 
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1 1 2 

Figure 2.1: for a particle in an infinite square well, such that Tr = SOTd. The distribution of 

coefficients is Gaussian, with n = 40 and an = 2. For the form of the square well that I am using, 
see Section [A.l.l[ 

1 1 2 

Figure 2.2: ||^|| for a simple harmonic oscillator. The distribution of coefficients is Gaussian, with 
n = 6 and cr„ = 2. For the version of the simple harmonic oscillator that I am using, see Section 

kJ3. 



A 



/ 1 




Figure 2.3: Plots of for a particle in an infinite square well at t = 0, t = Tr/2, and t = Tr/4. 
Note that the scale is different for the rightmost plot - though conservation of probability shrinks 
our peaks, they have the same shape as the peaks in the other two plots. 



CHAPTER 2. QUANTUM REVIVALS, FULL AND FRACTIONAL 



21 



happen at rational fractions of the revival time. 



2.3 Fractional Revivals 

We have seen the existence of revivals - times when the original wavefunction perfectly or near- 
perfectly reassembles itself. We have also seen that in at least one case we can find revival-like 



phenomena at t <C Tr. Indeed, if we look at Figure |2^, we might suspect that something interesting 
happens at most rational fractions of Tr. It turns out that while perfect revival^ at times other 
than the revival time are rather unusual, an initial wavepacket will give rise to arrays of other 
packets (which are not, in general, identical to the original packet) at rational fractions of Tr. The 
motivation for this is more algebraic than physical - it has to do with our ability at rational fractions 
of Tr to group together states with the same phase. This phenomenon was treated quite definitively 



in [|Averbukh and Perelman, 1989|] , and the next section will closely follow their derivation. 



2.3.1 The General Case 

In this section, we will restrict ourselves to wavefunctions that can be well-approximated by the 
following: 

'^{x,t)= ^ ci^k+n)^{k+n){x)exp-27:iik— + k'^—j,k = n-n (2.15) 

fc=-oo \ d Rj 

This is a fairly large class of problems - virtually any problem with a well-localized distribution of 
Cn around some mean value n, and any problem in a potential with a quadratic spectrum. This 
derivation makes no assumptions about the dimensionality of the problem, but does assume that 
the spectrum is a function of only one quantum number. Though we don't believe the generalization 



of these exact fractional revivals, se e [ p^^oinaz and Newman, 1999 1 



See [Loinaz and Newman, 1999 1 for a discussion of such revivals. 
'^We apologize in advance for the lower limit on the sum, which should only go to —n. The idea is that while 
all the k < ~n terms should have weightings of exactly zero, we can make the approximation that their weighting 
coefficients are very small. For somee sums and integrals these limits will greatly simplify our calculation. We've also 
dropped the leading exp {—2-KiEnt) term, as it has no impact on the relative phase of the various eigenfunctions and 
thus no physical significance. 
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Figure 2.4: Fractional revivals in the infinite square well. Our initial wavefunction is a Gaussian 
packet, centered on x = L/2, with = 0.003L. The symmetry of our wavepacket allows it to have 
a full revival (modulo an overall phase factor) at i = Tr/S. At times below that, though, we should 
be able to count the fractional revivals. 
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to more variables to be particularly problematic, we have not considered it here. 

We will look at the wavefunction at rational fractions of the revival time, t = ^Tfi, where 
p, q £ Z and p and q are relatively prim^. At such times we might expect that the t/Tji terms in 
the time evolution will reduce to some set of rational numbers, and that many values of k will have 
the same phase contributions. We will begin by looking only at exact times, but will find that our 
result is also a good approximation at t = ^Tr + At. 



If we now define 



mod 1, (2.16) 



which characterizes the deviation in phase of each eigenfunction from classical oscillation, we can 
write 

^ix,t=^TR)= CkMx)exp(^-2Tri(^k^^^+(l)k^y (2.17) 

Now, consider the conditions for some of these phases to match. If we can collect the terms of the 
sum into a handful of equivalence classes (defining two terms as equivalent when they have identical 
(pk), we can turn the infinite sum over k into a finite sum over another variable^. This is a question 
about the periodicity of (pk which can be easily answered. We are looking for the minimal I such 
that p,q G Z and for all fc. 



4>k = (t^k+i 

-fcVod 1 = -{k + lf modi. (2.18) 



This leads us to two conditions on /, 

^ G Z, (2.19) 

'-— G Z. (2.20) 
*This relatively prime business is Koing to fo rce the set of reviva ls that we can resolve to fit a Farey sequence. For 



more on this, see Appendix H or [Harter, 2001 1 or [Weisstein, 1999 1 



^Well, that finite sum is itself over another infinite sum, but since when did you get something for nothing? 
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There are two important solutions to Equation 2.19, I = q and / = q/2. Our task is to find out 



when the I = q/2 is consistent with 2.2C , and when the minimal value of I is just q. There are three 
cases here for us to consider, 



Case 1: q is odd 



For odd q, I = q. While I = q/2 would satisfy Equation 2.19| , q/2 is not an integer, as we require 



Case 2: q contains more than one power of 2 



When q contains more than one power of 2, Z = q/2. This obviously satisfies Equation |2.19| , and 



Equation 2.20 reduces to pq/4: G Z, which is also true. 



Case 3: q contains only one power of 2 



If q contains only one power of 2, / = q. While I = q/2 satisfies Equation 2.19| , it does not satisfy 
Equation 2.20| . That would require that pq/A G Z, but since p and q are relatively prime, p must 
be odd, and by assumption g/4 is half-integer. 



We have established that (pk is periodic, and we will be able to turn our sum over k into a sum 
with / terms. Note that if we were to consider times that were not rational fractions of Tr, we 
would be able to group the phase contributions from the k'^ term. With this relationship in hand, 
we may make a guess about how we may rewrite the wavefunction. Let us first define 

*d(^,i)= Ck^|Jk{x) exp ( -27:ik^Y (2.21) 

fc=-oo ^ ^^'^ 

a version of the wavefunction with the "dispersion" terms removed. We postulate that we can then 
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write the wavefunction as a weighted sum of time-shces[^ of ^ 



Cl J 



^(^x,t = ^Tr^ = J2 ^s^ci (^x, ^Tr + jT,i^ . (2.22) 
We will now demonstrate that Equation 2.22| is valid by explicit construction of the coefficients Ug. 



First, we plug Equation |2.21| into Equation |2.22| , 



^{x,t = ^TR) 



^Cs^'d [x,t+ -Tel 
s=0 

^ Cktpkix) exp f-2mk ^ is/l)Td 



=0 fc=— oo 
oo 



Cki^kix) exp ( -2mk— ) «s exp {-2TTijk^ . (2.23) 



We now compare the result with Equation 2.17 and arrive at the condition 



exp {—27ri(pk) = exp ^—2TTi jk 

s=0 



(2.24) 



for k = 0,1,..., I- 1. 



We now need to find an explicit expression for the constants a^. If we multiply Equation \2.24 
by exp (27rzyA;), where m is an integer, and sum over k, we can use the completeness relation of 
the Fourier sum to extract the coefficient: 



^"Recall that the rational numbers form a dense subset of the real numbers. This means that we can approximate 
the wavefunction arbitrarily well at any value of t with slices of ^d, which in turn implies that the set of time slices 
of '^ri taken at rational time s form a basis for the Hilbert space of a particular problem. This was first observed in 
[Aronstein and Stroud, 1997 1. 
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i-i i-i i-i 



fc=0 s=0 
l-l 

l^as6s,m, (2.25) 

s=0 

lgg-2^.(<^.-Tfc). (2.26) 



k=0 k=0 s=0 

l-l 



I 

k=0 



We have constructed the as, but we can rewrite equation ( 2.26| ) in a form more conducive to 



computation. First, we define how we will step through the possible values of m: 

vl 

m' = (m + 2— ) mod L (2.27) 



With this in place, we are prepared to derive the following: 



ttm' exp ( — 27ri ( ~ H — 



from which, shifting the summation index by 1, exploiting the periodicity of 4>k, and rearranging 
the equation for convenience, we conclude 

am' = ame y^^"). (2.29) 

It is easily confirmed that when q is odd or contains more than one power of 2, the a^'s with 
m = I and m = 2^ will have opposite parity, and all of the will have the same modulus. When 



with even and odd indices will have different moduli. As it happens, the coefficients with even 



q contains only one power of two, both for m = / and m = 2^ will be even, and the coefficients 
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indices are all zero. 

Now, let us compare our result with the actual wavefunction at t = ^Tr + At. The actual 
wavefunction is 



^ (x,^TR + At) = V CkMx)exp-2Tri(^'^k + 



—k + -k + —k 
Tci q Tr 



Using Equations 2.23 and 2.24, we find that 



(2.30) 



^' ( X, t = -Tr + At 



Y,as'^d(x,^TR + At+jT^) (2.31) 

Qo l-l / s \ ( Tr At \ 

^ ^ exp y-^m-k^ CfcVfc(^) exp -27ri jT"^ + jT^^j[2.32) 

(2.33) 



fc=— oo s=0 

oo 



Cktl^k{x)exp-2TTi f-'^k + ^k + -k'^ 
J-^ \qTd Tci q 

k=—oo 



which will be a good approximation so long as kAt/TR <^ At/T^i, where we must consider the largest 
"relevant" k. This is, perhaps, a closer correspondence than we should really expect - remember that 
at the start of this derivation we only required that our rewriting of the wavefunction be equivalent 
to it at t = [p/q) Tr. 

We may ask what this means for || 'I'll . At a rational fraction of the revival time, we have 



^ ( x,t = ^TR 

q 



J2 asal^ci (x, ^Tr + jTa) Ki fx, -Tr + jT,i 

s,z=o ^ ^ I- / \ q ' 



(2.34) 



a product of several time-slices of "ifd- In cases where we didn't have a well-defined wavepacket 
to begin with, or in the case where ^'cz doesn't preserve this wavepacket, we should not expect 
discernible fractional revivals. However, to the extent that each slice of ^'c^ contains a peak that 
has no significant overlap with any of the other time-slices in question, the cross-terms in the sum 
will make a negligible contribution and we will have several well-defined wavepackets - a fractional 



revival. To see this, consult Figure 2.4, where you should be able to count the fractional revivals. 
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Of course, we have (as of yet) no guarantee that is really free of dispersion, or that wavepack- 
ets really follow their classical paths. Though we have attached the label of "classical" to ^'c/ and 



Tci, we have yet to prove anything about them - that will have to wait for Section 5.5, when we've 
developed a technique for analyzing interference in quantum carpets. Before then, though, we can 



present some highly suggestive pictures, Figures ^ and |2.6 
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Figure 2.5: V{ar) (above) and (below) for a particle in the Poshcl- Teller potential (see 

Appendix ^). The dotted line in the above plot indicates En, in the lower plot they indicate four 
different classical paths of energy En. The Poschl- Teller potential has a purely quadratic spectrum. 
For more details on it, see Section A.2.3| . 



Figure 2.6: V{ax) (above) and ||^c/|P (below) for a particle in the Morse potential (see Appendix 
^). The dotted line in the above plot indicates En, in the lower plot they indicate four different 
classical paths of energy En. The Morse potential has a purely quadratic spectrum. For more details 



on it, see Section A. 2.1 
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2.3.2 A Closer Look at \l/cz and 

If you take a moment to write down the particular am for some p/q, you will find that all values in 
question appear twice|^. For example, at p/q = 1/5 we find 



ao = ao, (2.35) 

a2 = ao exp ^— 27ri-^ , (2.36) 

a4 = ao exp ^— 27ri-^ , (2.37) 

ai = ao exp ^— 27ri-^ , (2.38) 

as = ao exp f — 27ri- j . (2.39) 



Of course, we can find a reason for this. Let us suppose that from Equation p.27| we can construct a 
function m(n) such that m(0) = and m(— n) = m{l — n). If this is true, we can further simplify our 
expression for fractional revivals. Obviously, m(n) would depend on whether I was odd, a multiple 
of more than one power of two, or an odd multiple of two, so we will have to perform this analysis 
for three different cases. Fortunately, the proof that these a^'s come in pairs follows the same 
outline in each case. First, we will confirm that (m(n) +m(— n)) mod / = 0, which will tell us 



that afn(n) and Om(-n) connect ^'c/ {x,t + At) and ^ci — At) (see Equation 2.22 ). Second, we 
will prove that am(n) = o,m(-n)i to see if we could further simplify the sum in Equation 2.22 . This 
is to be done by induction - since a„(o) = a^i-o), if am(n) = «m(-n) implies am{n+i) = am(-n-i) 
then am(n) = am{-n) for all n. 

From Equation 2.2g| , we know that 



— h - j , (2.40) 



While the explicit examples that they calculated reflected the results of this section, Averbukh and Perelman did 
not state or use them. 
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and if we change the sign of re and shift the indices by one we can rearrange this to read 



' m{—n — 1) p 
am(_„_i) = a^(_„) exp -2m [ ^ h - 



(2.41) 



Using our inductive hypothesis, am(n) = 0'm{-n)) we need only show that 



-m(—n — 1) m(re) p 



(2.42) 



If all of this holds, we can define 



2 , / odd, 
L I even. 



and rewrite Equation ^.22 as 



(2.43) 



(2.44) 



Let's examine the three cases: 



Case 1: q is odd 



If q is odd, then I = q, and m' = (rre + 2p) mod q, from which we can deduce 



m(re) = 2pn mod q, 

(^m{n) 



n ■ f P , 2pre mod q 
ao exp ivrz n — h >^ 



(2.45) 
(2.46) 



The definition m(— re) = —2pn mod g is consistent with our condition that m{—n) = m{l — re). 
Moreover, under this definition rre(re) + m{—n) = 0, satisfying our first condition. Turning our 
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attention to Equation 2.42, we see that 



-m{—n — 1) 
I 



-{—2pn — 2p) mod q 



2pn 2p 

mod q -\ 

Q q 

I q 



(2.47) 



and our theorem holds. 



Case 2: q contains more than one power of 2 



If q contains more than one power of 2, then / = q/2, and m' = {m + p) mod q/2. We may then 
say that 



m(n) 



pn mod q/2, 



n-l 



o . / P , sr^ pn mod q/2 \ 
«m(n) = oo exp 2Tn \n-+2_^ — I . 

\ ^ N=0 ) 



(2.48) 
(2.49) 



The definition m{—n) = —pn mod q/2 is consistent with our condition that m{—n) = m{l — n). 
Moreover, under this definition m(n) + m{—n) = 0, satisfying our first condition. Turning our 



attention to Equation 2.42, we see that 



-m{—n — 1) 



— {—pn — p) mod q 

pn , 2p 

— mod g H 

q/2 q 

min) p 



(2.50) 



and our theorem holds again. 
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Case 3: q contains only one power of 2 

If q contains only one power of two, then again I = p. The principle difference from the (7-0(1(1 case 
is that uq = 0, while ai is non-zero. The proof in the g-odd case holds for the odd m, while for even 



m, ttm = 0, and we may easily define a;_ 



0. Thus our theorem holds in this final case. 



2.3.3 The Infinite Square Well 



Equation p.44| is particularly useful when we can find some connection between the paired slices 
of "ifd- The eigenfunctions of the infinite square well are harmonic functions, particularly easy to 
combine with the time evolution term in both and "^ci- Perhaps the first thing to do is plot a 



particular example - Figure ^ shows what looks like two wave packets bouncing off the walls of 
the well. 




t-> 



Figure 2.7: "^ci for a roughly Gaussian spatial distribution, centered at rc = a/3. 



We start out from equation (|2.21| ) and fill in the specifics of the problem: 



'^d{x,t)= ^ C(fc+fi)W--sin 



fc=— oo 



(A; + n) TT 



X exp —2Trik— 
L J \ Tci 



(2.51) 
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We then rewrite the sine terms as a sum of exponentials, 

>^ [2i { kir f 2L \ kn f 2L \ \ 

fc=— oo 

= E ^^^^12 (^i"x(^-i*)+sm^(x + §t) (2.52) 
+i cos ^ (x - ft) - z cos ^ (x + ft 



Note that f is the velocity associated with classical motion in the square well. If we plug in 
t = lb jTci, we find that we can rewrite the arguments of the trigonometric functions as 



_^ 2^Th A 2Lf . Now, inserting Equation |2.52| into Equation WM, we find 

Q -'- cl / >> * ' ' ' 



Prr. ^ [2 f kTT f 2Ls^\ kTT f 2Ls^ 

h 

= ^as(^'(x- Axs,0) + + Axs,0)), (2.53) 

s=0 

where we have defined Axg = 2Lj. 

It is rather surprising that in the infinite square well, at any rational fraction of the revival 
time, the wavefunction can be treated as a weighted sum of translations of the initial wavefunction. 
Note that \I'(— x,0) = ^'(2L — x,0) = — ^'(x,0), which ensures that each pair of translations will 
meet the boundary conditions of the square well - indeed, these translations can be thought of as 
disturbances in a dispersionless string, with the leftward (rightward) translation contributing the 
vr-phase-shifted reflection of the rightward (leftward) translation. 

The result that all fractional revivals in the square well consist of translated copies of the original 
wavefunction was stated but not proven in [Aronstein and Stroud, 1997|. There, the analogy with 



the dispersionless string is simply asserted, and the translations derived from that. 

This impressively simple result is, of course, a direct consequence of the overwhelming simplicity 
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of the infinite square well. That does not mean, however, that we cannot gain important physical 
insight from this result. We can think of "ifd as not changing shape during its time-evolution because 
there are no features on the bottom of the well. We might then guess that the degree of dispersion, 
and thus the degree to which the fractional revivals do not resemble the original wavepacket, depends 
on the features of the potential. 

2.3.4 A Final Note on "^ci 

Obviously, "ifd does not satisfy the Schrodinger equation - can we find a similar equation that it 
does satisfy? We can examine a few derivatives, and see if we can reconstruct something analogous 



to the Schrodinger equation. First, referring to equation 2.21 , defining E = |£'^|we can calculate 
the following derivatives: 



5^ 



- j;cfeVfc(-i^fc)e-^^'=* (2.54) 



dt 



Invoking the Time-Independent Schrodinger Equation, Ekipk = ~^g^ + ^^fci we can find 



k 

= V-^d - ^ CkEkipkC 



Ekt 



k 



V^ci - -Y^Ck {Ek - Ek) ^ke-''"'\ (2.56) 



dt 

k 

which we can rearrange to the following Schrodinger-like equation: 



^ = + V^d -Y.^k {Ek - Ek) ^,e-^^^ (2.57) 
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Brief examination of Equation 2.57 should reveal why this did not prove to be a successful vein of 
analysis. We have lost linearity, one of the most attractive features of the Schrodinger equation. 



2.4 Summary 

One of the principal differences between classical and quantum mechanics is that quantum systems 
have more than one time scale. Because of this, we can find many interesting phenomena in their 
time-evolution. In this chapter, we saw that the phenomenon of full and fractional revivals are quite 
general consequences of weighted sums of complex exponentials. We have also seen that we can 
express fractional revivals of a wavefunction as a weighted sum of a "classicized" wavefunction, ^ 
regardless of the original spectrum, we give ^'cz a linear spectrum, greatly simplifying its evolution. 
A new relationship between the weighting coefficients in that sum, am and a_m, was proved, leading 
to a proof that fractional revivals in the infinite square well are simply translations and refiections 
of the original wavefunction. 



Chapter 3 



Quantum Beats 



The object of consideration in this chapter^] is a weighted sum of exponentials, 

/(i) = ^Pne*^(")*. (3.1) 
n 

As one might expect, these sorts of sums exhibit revival phenomena just as in the previous chapter, 
when the Pn were functions. In the literature, this sort of sum is said to produce "quantum beats," 
and they are important because many laboratory measurements are of signals that can be written 
in this form - for example, the fluorescence of an ensemble of atoms, each excited to one of a few 
almost identitical energy levels. This is a useful enough phenomenon to have given rise to techniques 
such as quantum beat spectroscopy^. Because this is effectively a zero dimensional problem, there 
are fewer dynamical features for us to consider. Because of the weighted exponential form, we will 
see revivals and fractional revivals. What is interesting about this case is that we will not only note 
that they exist, but find a way to characterize how long the smooth revivals live before dephasing. 

As in the previous chapter, we are primarily interested in distributions of the weighting constants 
Pn centered around some mean value, n, with spread An, such that 1 <^ An <^ n. In that case we 



^Our treatment of quantum beats is based largely on that in [ Leichtle et al., 1996b , Leichtle et al., 1996a 
^For more on the experimental aspects of quantum beats, see [Silverman, 1995|. 
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can define k = n — perform a Taylor expansion of E{n) around n, and define 

2tt _ E^^\n) 



we can write f{t) as follows: 



(3.2) 



fit) = exp (lEnt) Pn+k exp 27ri (^^/c + ^/c^ + + ...^ . (3.3) 

Henceforth, we will ignore the leading phase term, as it has no impact on ||/(t)|p. Note that the Tj 
are not necessarily positive. When we interpret them as characteristic time scales, we will consider 
their absolute values, but allowing them to be negative simplifies our formalism. We will use this 
form of f{t) for the remainder of this chapter. For examples of two different sets of P„, one a 
Gaussian distribution and the other a "top hat" distribution, both with T2/T1 = 200, see Figures 



3.1 and . These specific examples will be treated in detail later in this chapter. 



3.1 The Early Phase of the Evolution 

The convenience introduced by considering P„ instead of ipnix) is that we can use the Poisson 
summation formula to recast the sum in Equation 3^ as a sum of time-separated signals. The 
Poisson summation formula]^ is 

^ / (7(A;)exp(-27riZfc)dA:= 9k, (3.4) 

I J —00 J 

where g{k) is an arbitrary continuous function such that g{k) = g/.. Using this formula, we can 



rewrite Equation as 
fit) 



00 

E 



I dkPik)exp27ri(^(^-^^-l^k + ^^k^ + ^^k^ + ..}j , (3.5) 



'For a derivation, see |Courant and Hilbert, 1953 | 



CHAPTER 3. QUANTUM BEATS 



tElL 




t T, 



100 150 200 




1 T, 



100 



fHL 



fHL 



1 

0.8 
0.6 
0.4 

0.2 



1 2 3 4 5 6' 



1 

0.8 
0.6 
0.4 
0.2 




fHL 



fHL 






100 



101 



102 



tTi 




Figure 3.1: A Gaussian distribution of the P^, with T2/T1 = 200 and cr„ 
early evolution, by Equation 3.29 dephasing should occur at t « 3.125Ti. 



8. In the plot of the 



Figure 3.2: A Gaussian distribution of the Pk, with T2/T1 = 200 and N = 8. In the plot of the 
early evolution, by Equation 3.53 dephasing should occur at t « 3.125ri. 
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where we interpret each integral term as a signal in time. Of course, this only constitutes a useful 
simplification when each integral term has a width in time that is less than the distance from its 
neighboring signals. Also note that our choice of P{k) has a significant impact on the tractability of 
the integral. Fortunately, the only restriction on P{k) imposed by the Poisson summation formula 
is that it be continuous and that P{k) = P^. 

3.2 Fractional Revivals 

It would be surprising if we were not able to combine our use of the Poisson summation formula 
in the previous section with the fractional revival technique of the previous chapter. As there, we 
might begin by assuming that it is the influence of T2 that will dominate most of the behavior 
beyond the scale of Ti. Our definition of a fractional revival must change a bit - here we mean the 
appearance of smooth oscillation at times less than T2. To that end, we will start near a time that 
is an integer multiple of Ti and also close to a rational fraction of T2, 

tp/g = ITi + At = ^T2 + ep/qTi + At. (3.6) 



Note that |ep/q| < 1/2. At t = tp/g we can rewrite Equation ^ as follows: 



f{t = tp/q) = Yl (^^^) + (^P/-? + ^) + f ^ + ^ ] ^k^ + • 



(3.7) 



If we define 



Wk = exp27ri ( -k'^ ] , (3.5 



Sk{At) = p^^i^exp27Ti(^^k+(^ep/g + ^^^^k^+(^l + ^^^^k^ + ..}j , (3.9) 
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we can rewrite the sum as 



f {t = tp/q) = ^ WkSk (At) . 



(3.10) 



fc=— oo 



Note that Wk, the weighting factor, is related to (j)k (Equation 2.16 ) from the previous chapter 
Wk = exp (— 27ri) cpk- In analogy with that, we define 



q, q contains one or zero powers of 2, 
|, q contains more than one power of 2, 



the minimum period for Wk, just as we defined I, the minimum period of (fik in Section Since 
we can write = Wk+mj, where m is an integer, we can use the formula 



oo J — 1 oo 

^ «A; = ^ ^ ar+mj 
k=—oo r=0 m=—oo 



(3.11) 



to simplify the sum in Equation |3.10| . The periodicity of Wk lets us extract the weighting factor, 
and rewrite Equation ^.10 



as 



j-l oo 
f {i = tp/q) = ^Wr ^ Sr+mj (At) . 



(3.12) 



r=0 m=— oo 



This is a form to which we can apply the Poisson summation formula. 



i-l oo „oo 

f {t = tp/q) = Wr / s (r + mj, At) exp {—2iTimk) dm, 

r=0 k=-oo-^~°° 



(3.13) 



where s (r + mj, At) is, again, a continuous extension of Sr+mj (At)^. This is still not so simple a 

form as we might like, so we introduce a change of variables, x = r + mj, and finally complete our 

''We just applied the Poisson summation formula to a sum of functions instead of a simple sum. Fortunately, we 
can think of each At as indexing a separate sum, and thus what we have done is a shorthand for applying the Poisson 
summation formula to a large (uncountably infinite) number of different sums. 



CHAPTER 3. QUANTUM BEATS 41 
serious manipulations, 

j-l oo 



f {t = tp/q) = -'^^'^r ex.p ^iri-k^ J s {x, At) exp ^—2TTi—kJ dx. (3-14) 
We may, of course, rearrange some of these terms in order to make this equation look like more of 



an improvement over Equation |3.7| . We may group all of the time-independent terms together into 
one new weighting coefficient, 

Wfc = - Vexp27ri fr^- +r^V (3.15) 



and a time-dependent term, 

Sk (At) = j dxP (n + x) exp Ini (^(^^ - x + (^e^/^ + ^x^ + (^l + ^x^ + ...^ , 

(3.16) 

and arrive at our final form for this sum, 

oo 

f{t = tp/q) = Yl WkSkiAt). (3.17) 

k=—oo 



The form of Sk (At) should evoke the integral in Equation 3^. In fact, it arises by making the 
substitutions 

t \ /At k\ 

t f At\ Ti 

t ( , At\ Ti 

Where the original solutions are centered on t/Ti = — /, the fractional revivals are centered on 
At/Ti = k/q, and the T2 and Tj terms are shifted relative to the Ti term, though, due to the 
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hierarchy |Ti| ^ IT2I ^ IT3I <C the shifting should be small. Because of this, the fractional 
revivals should look very much like the early evolution, further justifying our calling them "Fractional 
revivals." That one is able to arrive at the equations for the fractional revivals just by making 
substitutions, without doing any new integrals, is quite a computational nicety. Finally, note that 
though each constitutes a fractional revival, as in the previous section these fractional revivals 
are only distinguishable when they don't overlap significantly. 



3.3 A Special Case: The Gaussian Distribution 



The integral that will be used heavily in this subsection is 



-ax +bx 



vr ,,2 



Re {a} > 0. 



(3.21) 



We will define P{x) to be 



P{x) 



\/27rAn2 



exp 



{X — n) 
2An2 



(3.22) 



a Gaussian distribution of P„. In order to simplify the integrals, we will assume^ Tj = 00, for j > 3. 

°For a treatment of a Gaussian distribution of coefficients with non- infinite time scales Ti, T2, and Ts, see 
[ Leichtle et al., 1996b| . The solutions are in terms of Airy functions, and while the results are quite interesting, 
plenty of interesting results come from an case with just two time scales. 
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3.3.1 Early Evolution 



We are able to evaluate the integral in Equation 3.5, since it is a product of two Gaussian functions: 



1 

dk— - 

oo V27rAn2 



exp 



A;2 \ 



2An2 



1 exp 



27ri 



— - / k + —k' 
Ti T2 



dk exp I 27rz [ — l \ k 



1 



vr 



1 



■ exp 



r2 2An^ 

exp 



- 4TTiAnH/T2 \ l-4TTiAnH/T2\Ti 




(3.23) 



We can write this as a product of two Gaussians, one of real and one of imaginary argument, by 
rewriting the argument of the exponential, 



l-bi 
27r2An2 



+ 



abi 



1 + 62 1 + 62' 

27r2An2 



l-'i-KiAnH/T2 l + lQn^AnH^/T^ I + Wtt^AuH^ /T^' 



(3.24) 
(3.25) 



defining the (time-dependent) widths, 
and writing f{t) as 



1 



t 



2 1 



^ , ^ + 8An^^ 



1 



(3.26) 
(3.27) 



oo 

, A/l-47rz 



Z=— oo 



y^l -47riAn2t/T2 



exp 



2a2(t) 



exp 



2^f(t) 



(3.28) 



Equation 



3.28 



holds exactly so long as our assumptions are true (note that a Gaussian distribu- 
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tion of coefficients is always going to be truncated, so long as there is some minimum n, introducing 
an error function). f{t) consists of a sum of Gaussian packets separated in time and expanding as 
time increases, multiplied by a phase factor. These Gaussian packets begin to interfere significantly 
when their widths becomes comparable to their separation, 2 x 2ar{t) = 1, 



1 t2 

+ 8An2— = 1, 



V2 V 27r2An2 T^ 



leAn^ — 



t \^ 1 1 



T2J 4 TT^An^' 

' 1 



T2 8 An V vr^An^ 



t T2 1 1^ 1 



Ti Ti 8 An V 7r2An2 

^ (3-29) 



Ti 8 An V 27r2An2 

Notice that this dephasing time is roughly proportional to 1/An. If we then write At « T2/2An as 
the uncertainty of this state in time and consider 

AE = h{E{n + An) - E{n = An)) 

^ ( 2An 4nAn 

= 27r/i \ 

V Tx T2 

/I 2n\ 

= 47rnAn + — J , (3.30) 

as the uncertainty in energy, than we have recovered an uncertainty relation, 

AtAE ^ 2-Kh + 2n^ . (3.31) 

I believe this uncertainty relation to be related to one of the same form that appears in textbooks 
on quantum mechanics^. There, At is the time it takes for the expectation value of an observable to 
change by one standard deviation - a measure of how long it takes a system to change substantially. 



See, for example, [Griffiths, 1995, 112-114] 
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What we are measuring here is not stated in terms of observables, but it is a way of measuring how 
long it takes our signal to undergo an important, qualitative change. 

3.3.2 Fractional Revivals 



Fortunately, to consider the case of fractional revivals we need only apply Equations \i.l8 and 3.19 



to Equations |3.28| . Doing this, we find 



fit = t 



p/g) 



E 

k=—oo 



1 - 47riAn2 (ep/^ + At/Ti) Ti/Ta 



: exp 



At I k 
Ti ^ q 



2cr2 [t = tpjq) 



exp 



Af I k 
Tl ^ q 



^'^^Ui = tp/q) ] ' 

(3.32) 



(7, 



{t = tp/q) 



1 



(3.33) 
(3.34) 



Unfortunately, we cannot simply apply Equations S.lg and 3.19| to Equation |3.29| , as the separation 
between peaks is no longer 1, but 1/q. Instead, we require 2ar [t = tp/^j = 1/q. This leads us to 
the dephasing condition 



16An' I fp/, + — 1 



A A 
A A 

^p/q + ^ ) 



1 1 

T2 1 FT 



1 



Tl 8An V g2 7r2An2 ^ 
1 1 



-^ 1 



1 



Tl 8An g 
gTi 8 An 2g27r2An2 



q'^ir'^An'^ 
1 



An > 1(3.35) 



Comparing Equations 3.35 and 3.2£ , we see that for a fractional revival at t = tp/g, we will see 
almost exactly the same time evolution we saw at early times, contracted by a factor of 1/q. Refer 
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again to Figure 3.1 if you don't believe us. 



3.4 A Special Case: The Flat Distribution 

The integral that will be used heavily in this subsection is 



TV 

e~ 

-N 



the integral of a Gaussian truncated at itA^, where Erf (x) is the error function]], defined as 



2 / 



Erf (x) = —= e~' dt. (3.37) 

Erf is odd, Erf (0) = 0, Erf (cxd) = 1. Prom this we can conclude that as A'' becomes large, the error 
function term should go to one. We will define the "top hat" distribution P{x) to be 



P{x) 



1, |x — n| < N, 

(3.38) 

0, \x-n\ > N, 



which describes a flat distribution of width 2N centered around n. This is an interesting case for 
several reasons. First, it is the other obviously integrable case, along with the Gaussian case. It is 
also an example of a distribution with sharply limited extent in n-space, unlike the Gaussian. It is, 
in a sense, the most un-Gaussian of localized distributions. Finally, if we study limAr_^oo, we can 
gain some understanding of an "ultra-localized" spatial distribution]^ As in the previous example, in 
order to simplify the integrals, we will assume Tj = oofor j > 3. 



^For a discu ssion of Erf, see either the more traditional |Abramowitz and Stegun, 1965 ] or the hip and modern 
llWolfram, 200^ . 



^R,emember that our ultimate goal is to relate these results to quantum mechanics problems, and we would expect 
that a broad distribution in energy-space would correspond to a narrow distribution in position-space. While in 
doing this problem we will break a great number of mathematical rules. P{x) is not continuous, as the Poisson 
summation formula demands, which casts a shadow on any analysis that we do. Studying limjv^oo has us implicitly 
using negative energies, and if we were to try this in quantum mechanics we would have a non-normalizable state 
that probably wasn't differentiable, and thus not even a solution to the Schrodinger Equation. All the same, we can 
hope through this savagery to gain some physical insight. 
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3.4.1 Early Evolution 



As above, we are able to evaluate the integral in Equation 3.5, since it is a truncated Gaussian 
function: 



fit) 



N 



1 



N V2ttN'^ 



exp 



t 



I 1 27riA; + 2m — k 

T2 



1 



vr 



V -27Tit/T2 



exp 



87rit/T2 



2 r^^i 



-Erf 



2y/-2mt/T2 



(3.39) 



We can simplify the arguments of the error functions, 



-A-KijrN ± 2m (^^ - I 
2^J-2mt/T2 



t 



t VTi 



(3.40) 



and recognize the sum of error functions, which shape the wave packet for a given I, as being of the 
form 

Erf {a(j){A - x)) + Erf (a0(A + x)) , (3.41) 



where 



a 



A 



vr t 
2T^' 



2N, 
t \Ti 



(3.42) 
(3.43) 
(3.44) 
(3.45) 



From here, it is helpful to first examine Figures |3^ and |3^ , plots of Erf {(ftx) and Erf {—(f>x). These 
suggest that if we are to study the modulus of Equation 3.39| , that we would do well to consider the 



CHAPTER 3. QUANTUM BEATS 



48 




Figure 3.3: Re {Erf (^x)}, x G [-10, 10]. 
Figure 3.4: Im {Erf x G [-10,10]. 

real and imaginary parts separately. Doing this, we find 

\\Eif {a(t){A-x)) + Elf {a(t){A + x))f = ||Erf - a;))f + ||Erf + x))f + 

Erf {a(l){A - x)) Erf {a(f>* {A + x)) + 

Erf {a^* (A-x)) Erf {a(t>{A + x)) 

Re {Erf (q:^(A - x))}^ + Im {Erf (q!^(A - x))}^3H6) 
Re {Erf (a(/)(A + x))}^ + Im {Eri{a<j){A + x))}^ + 
2Re {Erf (a^(A - x))} Re {Erf (a0(^ + x))} + 
2Im {Erf (a0(A - x))} Im {Erf + x))} . 

Here we see that, roughly, Im {Eif {(f)x)} = + e(x) and Re{Erf = Sign(x) (1 + e'(x)), where 
e(x) and e' (x) are functions that describe the oscillations. Obviously, both e(x) and e'(x) decrease 
as jxj increases. We can then observe that in the region x < —A, the contribution from Re terms 
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= Re {Erf {a(l){A + x))}^ + Re {Erf {a(t){A - x))}^ + 2Re {Erf (a(/>(A - x))} Re {Erf (a(/>(A + x))} 
= (-l-e'(x))' + (-l)2 + 2(-l-e'(x)) 

« l + 2e'(x) + l-2-2e'(x) = 0. (3.47) 

The same calculation can be performed for the Im{} terms, and for x > A, thus demonstrating that 
anything interesting that happens in the modulus of this sum of error functions happens between 
—A and A. Although our characterization of the width of the pulse described by Equation 3.46 is 
not so transparent as the standard deviation of a Gaussian curve, we can say that we have described 
a pulse are centered on x = 0, with width 2A. 



Analysis of Equation 3.40| , the argument in our particular problem, is also not so simple as it 



was in the case of the Gaussian distribution. First, let us find the edges of our pulse - times that 



satisfy x = —A and x = A, respectively, in Equation |3.46| . We will define these two times, t2 and 
ti, as follows: 

^ ^ (3.48) 

^ = (3.49) 

Ti 1-2N^ 

Clearly, t2 — ti defines the width of the pulse, and (ti + t2)/2 defines its center. If we define the 
spacing between pulses as the distance between their centers, then that spacing is 

1 / 1 

Hl + 2iVft^l-2iV|J-,_^2iV|)'' 

We are now looking for the spacing for the the value of I for which the spacing between packets 
equals the width of a packet, the value of I for which we think the packets will overlap enough to 
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interfere substantially: 



^1 1-(2N^^ 
8NIT1/T2 1 



2 • 



I = (3.51) 



Ti 8iV 



We convert this into a time by multiplying it by Ti 



T2 



t = T,l = ^, (3.52) 



and finally write 

t _ 1 
T^ ~ 8N' 



(3.53) 



As in the Gaussian case, the time after which our signal dephases depends on Looking 
for an uncertainty relation, we have At = T2/2N and, as in the Gaussian example, AE = 
A-kKN (1/Ti + 2n/r2)and we recover an uncertainty relation, 

AtAE = 2'Kh + 2n^ , (3.54) 

the same relation that we found in the Gaussian example, though in that case the result was not 



exact. If we refer to Figure |3.2| , though, we find that our dephasing time doesn't seem to fall in the 
right place. I submit that the noise between the t = 2Ti and the t = 3Ti pulses is not the dephasing 
that we were looking for, as it is followed by the discernible top of the t = 3Ti pulse. 
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3.4.2 Fractional Revivals 

Applying Equations 3.1^ and 3.1£ to the results of the previous section, we find our signal to be 



47r^ 



fit- ip/g) - 2vh^\/ -2ni{e^^^+At/n)in/T2) "^""P [ 87ri(.p/,+AVr')(Ti/T2) J ^ ^^'^^^ 

I Erf 



Erf 



-47ri(ej,/,+At/Ti)(Ti/T2)Af-27rj(^-|^ 

2^ ~2ni{ep/^+At/T^){Ti/T2) 
' At I 



-47rj(ep/,+At/Ti)(Ti/T2)Ar+27ri('^*-'= 



2y-27ri(ep/,+At/Ti)(Ti/T2) 



which could be simplified. More interesting is, of course, our dephasing condition. The method that 
we use to find this is the same one that we used in the previous section. 

We started with our simplification of the arguments of the error functions, in Equation B.40| , 
which we will convert, 



1 I T/ Tift \\ 1 I f AA /TA / #-t 



2V T2\ ' t \Ti J J 2^ y^''^ TiJ\T2j\ (ep/. + At/Ti) (Ti/Ta) 

(3.56) 

We now look for the roots of this, 



±2iV, 



1 /At V 

(ep/, + At/Ti) (T1/T2) " ~q, 

At k _ ^ oAr^l f , 



J 1 / L^L \ 



At / Ti\ 



k Ti 
- ± 2N—e 



T2 
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and label those roots, 



Ti 1 - 2iV§ 

Ti 1 + 2iV§ 



We're now interested in the width of the pulses, 



Ail AtA _ 47V4^ep/, + 4iV4^| 



l-f2iV|) 



2 



and the spacing between pulses. 



We now find when these are equal. 



(3.58) 
(3.59) 



1 - i2N^ 



+ . . ^..T, = : 79- (3-61) 



2fe I ^\ = 1 1 

- = — — --e/ 3.62 
q Ti 8N q ^'^ 

and see that for fractional revivals as well, the dephasing time is proportional to 1/N , and decreases 
as q increases. Again, we see that the behavior near t = will appear again, squished by a factor 
of l/q. In fact, this result corresponds as closely to the one in the Gaussian case (Equation |3.35 ) 



as our results for the early evolution did. It would appear that the principle difference between 
the signals generated by these two dissimilar distributions is that the sharp edges of the top-hat 
distribution create the same kind of high-frequency "ringing" that we would expect from a truncated 
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3.5 Summary 

In this chapter we studied a zero-dimensional system that exhibits "quantum beats," which is both 
of experimental interest, and, due to its low dimensionality, is a nice system on which to introduce 
our Poisson summation formula technique. The principal advantage of this technique is that it 
allows us to quantify to dephasing of an initial wavepacket. We have seen that, for quantum beats, 
fractional revivals are related to the early evolution of the packet by the simple substitution of a 
few terms. In the specific examples that we have considered, a Gaussian distribution of weighting 
coefficients and a "top hat" distribution, we have discovered nearly-identical dephasing conditions. 
If An is the "spread" of the weighting coefficients, then the initial dephasing time goes as T2/8An, 
and the dephasing time for a fractional revival at t = tp/q goes as T2/8Anq, and the fractional 
revival has the same number of "clean" oscillations as there were oscillations near t = 0. That 
these distributions are so different from each other suggests that many well-localized distributions 
of coefficients should demonstrate similar behavior. 



Chapter 4 

The Connection Between Beats and 
Carpets 

One shortcoming of the previous chapter is that it describes a zero-dimensional problem which, 
while interesting, lacks the spatial features that we would like to consider. The main advantage of 
the quantum beats approach is that it allows us to not only identify pseudoclassical behavior at full 
and fractional revivals, but also to associate a particular lifetime with those revivals. In both of the 
examples that we considered we found that the lifetime of the revivals was proportional to 1/An, the 
spread of the wavepacket in energy-space, but two cases hardly exhausts all of the possible localized 
wavepackets. Although we could think of revivals of a spatial wavefunction as revivals of each point, 
we know that even a simple distribution of weighting coefficients will become complicated when it 
is combined with the eigenfunctions evaluated at a particular point. How, then, can we connect 
this approach with quantum carpets? 
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4.1 Abuse of the Poisson Summation Formula 
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So long as we can come up with a continuous extension of a set of functions fm{x), we have a 
generalization of the Poisson summation formula, 

OO OO „oo 

E /™(^)= E / fim,x)e-^''''^'dm. (4.1) 

m=-oo l=-oo °° 

Not all sets of functions will have an obvious continuous extension, but we can find extensions in 
cases where our eigenfunctions are continuous functions of m and x. So long as we are able to make 
this extension, it is easy to apply the formulae derived in the previous parts of this chapter to an 
n-dimensional problem, by making the substitution Pk — > Ckipkix). This is our crucial observation 
in this section, but its merit will not be clear until we have examined at least one case. 



4.2 A Special Case: Gaussian Weighting Coefficients in the Infinite 
Square Well 



The attentive reader will have noticed that the integral hiding in Equation |4JJ is a bit intimidating. 
In order to be sure that it's even worth thinking about, we'll do it for the simpliest case I can think 
of - the infinite square well. We'll take advantage of two characteristics of solutions to the square 
well: the eigenfunctions are just sine waves, easily converted into exponentials, and the spectrum is 
quadratic in the quantum number, meaning that a Taylor expansion of the spectrum is exact and 
gives just two time-scales. All of our integrals will reduce to Gaussian integrals, which we can do. 

4.2.1 Algebra 

For this case, we will write our wavefunction as 
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and defined 



" (4.4) 



27r^n/L^ vrn 
27r 2L2 



7r2/L2 TT 

'El = 1. 

T2 2n' 



(4.5) 
(4.6) 



What follows is a godawful bit of algebra, and in order to make it easier to follow (and ensure that 
I get it right), I'm going to temporarily replace all of these compound constants with simpler ones. 
I'm also going to ignore multiplication by leading constants - don't worry, I'll put it all back in 
later. If you're not interested in this, you won't have missed much by skipping directly to Equation 



4.27. We will define 



^. (4.8) 
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and use these, along with Equation 4.6 to define the following 



a = 


1 


(4.9) 


2An2' 


b = 


nvr 1 T2 

— X = vrf , 

L 2Ti ' 


(4.10) 


c = 


-X = 


(4.11) 


d = 


t 

27r— = 27rr, 


(4.12) 


f = 


t Ti 
T2 T2 


(4.13) 


9 = 


27r/. 


(4.14) 



Note that these are all real quantities. Having defined these, we can rewrite our wavefunction as 

^(x,t) = y^^ exp [-ak"^) sin {b + ck) exp - i [dk + fk'^) . (4.15) 
fc=— 00 

We now rewrite the sine term as a sum of exponentials, 

^'(3;,^) = — Y^sxp (— a/c^) (exp z (6 + c/c) — exp — i (6 + cA;)) exp i (— d/c — /A:^) , 
fc=— 00 



(4.16) 



combine the exponentials. 



^{x,t) = ^^Jy X] {ew {bi + iic - d)k - {a + if)k'^) - exp [-bi - i{c + d)k - {a + if)k'^)) , 
fc=— 00 

(4.17) 

apply the Poisson summation formula. 



1 Fl [a °° /"°° 
" 2i V IV ^ y (exp (6z + i(c - d - 5)A: - (a + i/)A;2) - 

exp [-bi - i{c + d + g)k - (a + if )k'^)) , 



(4.18) 



/=— 00 ' 
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and do the integration, using the formula dkexp (— aa;^ + = a/vt/q exp (/3^/4a): 

1^ 4(a+*/) J J ■ 

We have, in a sense, solved the problem now, but we have gained little unless we can understand 
this result in a simpler way than we understood the original expression. So long as the terms of the 
sum are localized in their own regions of spacetimeQ, and the region defined by one value of / does 
not overlap significantly with the region defined by another value of /, we can make the substantial 
simplification of studying just one term in the sum to understand a particular region. There are, 
however, some algebraic simplifications that we will make first. 

We perform the square in the exponentials, 

^ i/f t ^ (»P - - (4.0) 

exp [~bi - ^ ^ j j . 

We then turn the difference of exponentials into a product of a sine function and another exponential, 

') = yi/ E - c^p^) f * + if^) ■ (4.21) 

V LV vr^£-^ V a + z/ V 4(a + z/) J \ ^i{a + if)) 

From here, we would do well to write the arguments of the exponential and the sine in the form 
x + iy, 

- f,f-t - <« - ^/)) - - |$T^ (-/ - »)) ■ 

(4.22) 

Now, we will consider making the assumption that we can ignore interference between 

^This isn't obvious, but it is something we can hope for, as the Poisson summation formula works kind of like a 
Fourier transform. 
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terms in the sum. While this is not obviously true, this whole business is worthless if it isn't. Don't 
worry, we'll check this assumption at the end, just to be sure. One of the advantages of doing this 
is that we can use the formula ||sin(x + iy)||^ = (1/2) (cosh (2i/) — cos (2a;)). Doing all of this, we find 



-oo 



(4.23) 

Looking at this, we may notice that the adding the argument of the cosh to the argument of the 
exponential would complete the square. Fortunately, we may write the cosh in terms of exponentials, 
giving us 

oo 

ll*(-,Of = -y|=^ E ((e-P-(a^i;Sj^f)+exp-(a(5S^))x (4.24) 



cosh ( a 



This is the simplest form of this equation that I have found. Noticing the ubiquity of + /^, we 
define 

then fill in our various constants to arrive at the following. 



V '=°^''(2A^ ) ) ) 

^ 00 

(1 _ eos. (gf - 4.^ar$±fl) sad, • 
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4.2.2 Interpretation 



We are now prepared to discuss Equation 4.27] . Most of the structure is provided by the two 



travelhng Gaussians, which both originate at ^ = 0, r = — /, then separate and travel along straight 
lines to ^ = 1, T = — Z lb 1/2. This both provides a pleasing picture of a packet completing one 
classical oscillation in one classical period (Ti), and tells us that if the packets in two neighboring 
terms have cr(r) < 1/2, those terms will have negligible interference^. All of this is modulated by 
the final, oscillating term, which provides "interference without cross terms." Note that because 
|cos| < 1 and sech < 1, the oscillating term will always fall between zero and two, and makes 
a secondary contribution. We have left things in terms of the ratio T1/T2 as a reminder that a 
principal difference between classical and quantum dynamics is that a classical system would have 
only one time scale (Ti), while a quantum system may have more - in this case, two. Note that by 



Equation the limit as n — > 00 is identical to the limit T1/T2 — > 0, and that this fixes the widths 
of the travelling Gaussians above. That is, in the classical limit our packet does not spread. 

When we are not in the classical limit, three things prevent our travelling Gaussians from 
behaving classically. First, there is a damping term in front of the Gaussians, though this isn't 
difficult to interpret - it enforces conservation of probability. Second, their widths do increase, and 
when those widths become larger than the width of the well we would expect significant interference 
between, for example, the Ith term and the (/ + l)th term. Finally, there is the oscillatory term. 
In the cosine term, the (T2/Ti)^ term guarantees rapid oscillations, and the other term serves to 
curve the paths along which the cos x sech term is constant. The oscillation of the cosine term is 
modulated by the slowly-varying sech term. Although there is no conventional charateristic width 
defined for this function, as there is for a Gaussian distribution, it is clear that when the argument 
of the sech function becomes large, our Gaussians will look fairly smooth, and the argument of sech 
is on the order of one, the unsuppressed cosine will produce oscillatory interference-like effects. 



We can see all of this by plotting several terms from the sum in Equation 4.27 and comparing 



it to a similar region of the full expansion of Equation (see Figure [O] ). One thing that you 
^Of course, the "no overlap" assumption is going to break down in some interval near C — 1- 
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won't notice unless you plot some of these yourself is that where Equation 1.27 is valid, it is more 
computationally efficient - there are no cross terms to calculate. For an example of this, see Figure 



4.2 



^ 2 




^ 2 




Figure 4.1: These density plots show a solution of the infinite square well in the region r G [0, 4], ^ G 



[0, 1] for An = 5, n = 30. On the left are the five relevant terms from Equation 4.27, on the right 



are the middle 10 terms of the sum in Equation 1.2 



4.2.3 A Quick Look at Sech 



One function that appears in Equation 1.27 and may be unfamiliar is sech, the hyperbolic secant. 



Defined as sech(y) = 1/ cosh(y), it has a peak at y = 0, sech(y) = 1, and limij^i^o^^ = 0. In fact, it 



looks rather like a Gaussian wavepacket. How much? Consider Figure |4.3| , which suggests that near 
y = the functions can be given similar characteristic widths. For example, choosing a constant 
a = ArcSech(l/e^) 2.69 gives us exp (— y^/2cr^) w sech{ay/2a), which becomes exact at y = 
±2cr. We can generalize this result still further. Suppose we want exp (— y^/2cr^) = sech {ay/Aa)a.t 
y = ±Aa. We will first find that 

exp(-^)=exp(-^V2). (4.28) 
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6.188s 
207.891s 
106.047s 



Figure 4.2: All of these show the region r G [0.5, 1.5], ^ G [0, 1] for An = 5, n = 30, and the time 
it took Mathematica to plot them. On the left is the I = —1 term from Equation 4.27| , the middle 
figure contains the middle 20 terms of the sum in Equation |4.2| , and at right is a plot of the middle 
10 terms from that same Equation. 
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Figure 4.3: A plot of exp ^— (dashed) and sech (a^) (solid), —ha < y < 5a. The constant 
a = arcsech (^) 2.69 is chosen so that the two functions are equal at y = ±20". 
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Comparing this with 

sech (^~^^ ~ sech(a) 

= exp {-A^/2) , (4.29) 

we can see that a{A) must be 

a{A) = arcsech (^e-^'/^) . (4.30) 

Note that the domain of arcsech is (0,1], the same as the range of exp (— a^^), so 0({A) is defined for 
all real values of A. 

Now, our final step of generalization. Suppose we want to find the condition for 

to hold at y = ±1^0"^ = ±2^o'b. This equation is clearly satisfied when 

(y-a) = ±iAaa, (4.32) 

(y-b) = ±2Aab, (4.33) 

where I have used subscripts on the ± symbols to indicate their independence. It would be un- 
fortunate if our goal was just to find the intersections of two given functions. That is, if we were 
searching for the solutions of 



<76 



exp ( - ) = sech ( ^ ) ' (4-34) 



where we were unwilling to vary any of the parameters. This will quickly lead us to the transcen- 
dental equation 

Aaa + o = a{A)ab + 6, (4.35) 
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which I don't care to solve. However, if we fix A and choose to solve for one of the other variables, we 
will have an ordinary algebraic equation with perhaps one unsightly transcendental number floating 
around. Since we are only concerned with equality at a point, we can even allow a, b, aa, andcjf, to 



be functions. In that case, the result in Equations [4.32| and |4.33| are quite useful. They let us figure 
out what parameter values to use to get these functions to intersect at, say, one standard deviation 
from their centers, which is exactly the what we need to do in order to find dephasing times. 



4.2.4 Dephasing of the Wavepacket 

If we recall the previous chapter, the advantage of this sort of analysis of the wavefunction is 
that it allows us to quantitatively understand when the wavefunction makes the transition from 
semiclassical motion to characteristically quantum mechanical motion. We could assume that this 
happens when the sech term overlaps the Gaussian closer to r = - then we need only decide 
what constitutes substantial overlap. Let us suppose that it occurs when some particular number of 
"standard deviations" of the sech function reach a particular number of standard deviations of the 
travelling Gaussian (we'll not specify on which side of its center), so that the functions are equal at 



that point. This condition may be more clear when written like Equations 4.32 and 4.33: 



T+ 1 + 



2 

(r + l) 



±^o-exp('r), 
-^0-sech('^), 



(4.36) 
(4.37) 



where we have defined (Tsech(''") and o"exp(T") in terms of cr(r), from Equation 4.25 



<^sech{r) 
0"exp(T) 



a{A) f An 
~A^ 
An 



IT 



Mr) 



IT 



1 \2 / Ti 



(4.38) 
(4.39) 
(4.40) 
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Note that the ± in Equation 4.36 allows us to consider intersection on either side of the center of 



the Gaussian, while the fixed sign in Equation |4.37| indicates our choice to consider only the region 
between t = —I and r = 0. Keep in mind that what we are doing is, in spirit, exactly what we did 
in the previous chapter. We have used the Poisson summation formula to turn a sum of functions 
over all space into a sum of localized functions, figured out how to characterize their widths, and 
now we are preparing to find out when those widths overlap, in order to discover when interference 
effects become really important to the time-evolution - when our wavepacket has dephased. We 



will now combine Equations [4.32| and [4.33| to eliminate r + I and solve for o"(t): 



- AUsechiT) + - = ±Aaexp{T), 

a{A) f An . .V £, , . An , , 

a{A) /An , .V ^ /An , A ^ 



(4.41) 



We may now use the quadratic equation to solve for (Anair)) /vr. 



An 



±lA±n 



vr 



(4.42) 



Since a{A) and A are always positive, we know that the square root term will always be greater 
than one. We now square this equation, to recover something in terms of o"^(r). 



— ) a (r) 



a{A) 
~A^ 



(4.43) 



which we can further simplify 



(4.44) 
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before inserting Equation 4.25, 



1 \2 / Ti ^2 



2An 



T2 



2 VAnaC^);* I + A^ S ^ A^ 



]_T2 

271 Ti 



V2 T2 



47rAn Ti 



1 / ttA^ 

2 I Ana(A) 



1 + 



a{A) 



± a/1 + 2 



^4^ 



ttA^ 



2An' 



^.45) 



Note that the r at which dephasing occurs depends (to first order) on 1/An, as in all of our other 
examples. To check that this works, let's compare the lower of the A = 1 terms against plots 



of Equations [4.27| and in Figure |4^. Note that the plots agree quite nicely in the left half 
4 



^ 2 




^ 2 




Figure 4.4: These density plots show a solution of the infinite square well in the region r G [0, 4], ^ £ 
[0, 1] for An = 5, n = 30, and the white line indicates the lower of the A = 1 solutions of Equation 
[4.45| . On the left are the five relevant terms from Equation 4.27 , on the right are the middle 10 
terms of the sum in Equation 4^. Note that the ridges in the left plot only appear above the line, 
and in the right plot that this holds for ^ < 1/2 - for ^ > 1/2, the interference between terms in 



Equation 4.27 becomes non-negligible. 



of the well, but that where different terms in Equation 



4.27 



overlap, i.e. near the right wall, our 
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approximation does not reproduce the appropriate interference effects. Nevertheless, Equation 4.45 
does predict where ripphng begins in the left side of the well. If we wanted to set different standards 
for overlap, we would just need to change the value of A. We could also treat fractional revivals in 



the square well, but the reader should be able to combine the results of Section and those of 
this section to quickly find and study them. 

4.3 Solutions in More Complicated Potentials 



Obviously, we have little interest in actually doing the integral in Equation |4^. At the same time, 
there are many interesting potentials that we might consider which are not the infinite square well 
- which, indeed, do not even admit solutions in pure exponentials. Our analytical approach to 
these rogues will be to use the WKB approximation to render much of the information about the 
wavefunction in convenient exponential form. 

4.3.1 The WKB Approximation 

The WKB approximation^ is a formalization of the observation that if a wavefunction has sufficient 
energy that it undergoes rapid oscillations in comparison to the "features" of the potential below 
it, it looks much like a free particle. It is useful in the same sorts of semiclassical problems where 
we would expect to find interesting revival and carpet phenomena. Without going into too many 
details, we define 



Pnix) = ^/En-V{x), (4.46) 
which will behave rather like a wave number. It is subject to the Bohr-Sommerfeld condition, 



Pn{x) = mr, (4.47) 



For a more complete treatment of the WKB approximation, see iGriffiths, 19951 , 274-292] 
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which enforces quantization of energy, xi and Xr are, respectively, the left and right classical turning 
points for a particle of energy En- We are then able to write the wavefunction as 

i^nix) ^ exp (i I pk [x') dx'] + exp ( -i [ pk {x') dx'] , (4.48) 

VPn{x) \ J J VPn{x) \ J J 



with the understanding that both p {x') dx' and 1 / y^p„ (x) vary slowly in comparison to exp (i p [x') dx') . 
C and D are constants of integration. This approximation of the wavefunction breaks down when 
En is not much greater than V{x), near the classical turning points. 

4.3.2 The WKB Approximation and the Poisson Summation Formula 

The WKB approximation has not quite solved our problem, having both left that unattractive 
^/VPk (x) term in from of our exponentials. We are, however, en route to a solution, 

^ roc 

^ix,t) = XI / Ckipkix) exp {-iEkt) 

= Y dk i^^e-Kp-ii- j''pk{x')dx' + Ekt + 2TTkl)+ (4.49) 

t=— oo \ 

exp -i {p Pk {x') dx' + Ekt + 2M) , 
y/Pk(x) 'J 

but the remaining integral is still less than inviting. Of course, judicious choice of the c/j and d^ 
will make matters simpler, but we can do more than that. We now perform a Taylor expansion 
of p Pk {x')dx', Ek, c(k) / \/ pkix) , and d{k) / \/ p}.{x) around A; = 0. For the terms that appear in 
the exponential, we would keep up to order k'^ - k^ if we felt like using Airy functions. Keep as 



many terms as you Uke in the expansion of 1/ ^Jpk{x). The resulting mess will be some Gaussian 
integral which will be easy to evaluate in closed form. We should find some sum of exponentials, 
trigonometric, and hyperbolic functions, including some which we can identify as containing much 
of the interference information, and find the condition for wavepacket dephasing. 

Our solution relies on three successive approximations - the WKB approximation, the applica- 
tion of the Poisson summation formula (which involves an extension of the c„, dn beyond k = —n 
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to fc = — oo and eventual neglect of cross terms), and a Taylor expansion of the remaining functions 
- and obviously won't work in all cases. There are doubtless examples which can be integrated 
that neither bear resemblance to the square well nor require that we resort to these approximations. 
However, we have here a quantitative way of finding out the lifetime of a wavepacket in an arbitrary 
potential well. 

4.4 Summary 

In this chapter we demonstrated that the Poisson summation formula technique, developed on zero- 
dimensional quantum beats in Chapter |3|, can also be applied to problems with more dimensions. 
We demonstrated this by using the technique to solve the square well, arriving at a very elegant 
sum of successive travelling Gaussians multiplied by interference terms. Prom this, we were able 
to derive a dephasing condition, though instead of finding a simple time we found a function r (^). 
What's more, it only works in the left half of the well, due to interference that we had to neglect. 
Fortunately, we found a r oc (T2/T1) /4An component in the dephasing condition that is strongly 
reminiscent of our dephasing results from Chapter ||. The case of fractional revivals is left as an 
exercise for the reader. 



Chapter 5 



Intermode Traces and Quantum Carpets 



5.1 Multimode Interference 

The Multimode Interference technique]^ is a devilishly simple way of studying the interference pat- 
terns formed by a quantum system. It allows us to understand both the existence of canals and 
ridges in the carpet and the role that the spectrum and potential play in generating carpets. We 
will begin by rewriting the probability density as a sum of multimode terms (defined below) instead 
of eigenfunctions. Starting with a wavefunction, 

oo 

^{x,t) = Y,CnMx)e-'^"\ (5.1) 



n=l 



we define a multimode term as 



Unmix, t) = ^ {dnmtpn{x,t)ij^{x,t)eX.p-it{En - Em) + (5.2) 
dmntPni^, t)'tpmix, t) exp -it {E^ - En)) , 
dnm = CnCm = '^■mrt) (5-3) 



The first sections of this chapter are based on iKaplan et al., 200(]| |. I have tried to make their work more clear 
and rigorous. 
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and write the probability density as 

oo 

\\^{x,t)f= J2 f^nm{x,t). (5.4) 
n,m=l 

This is not, at first glance, a particularly good idea. We have simply grouped the terms in the 
sum differently and not achieved any obvious simplification. The advantage is that we may simply 
add multimode terms - they already contain all of the information about how the eigenfunctions 
interfere. If we can figure out a way to tease that information out of the multimode terms, we will 
have been successful. 



5.2 Characteristic Velocities 

As in the previous chapter, we will use the WKB approximation to study semiclassical cases (see 



Section 4.3.1). Our approximate eigenstates are then 



tpn{x) ~ 5^"' exp (i [ pn {x') dx'] + — £^exp ( -i [ pn {x') dx'] , (5.5) 

y'Pnix) \ J J \/Pn{x) \ J J 



where C„ and are complex constants of integration, and our approximate multimode terms will 
be weighted sums of eight intermode terms, which we define as 

Lnm.{x,t) = ]-dnm^= , \ , , exp ifciZ ( / {±2Pn {x) Pm {x')) dx + {En - Em) t] . (5.6) 
^ Pn{x) Pm{x) \J J 

Since we are interested in lines of constant phase, we differentiate the argument of the exponential 
and look for its roots, 

^±l(^y" {±2Pn{x')±^Pm{x'))dx' + {En-Em)t^ = 0, 

dx En - Em 



dt ±2Pn {x) ±3 Pm {x) ' 

AcJ En — Em 

""-"-^Ak = ^'^En-V{x)±2^Em-ViiV^ 



CHAPTER 5. INTERMODE TRACES AND QUANTUM CARPETS 
Each pair of quantum numbers, {n,m), gives rise to four velocities Vnm- 
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5.3 Characterization of the Velocities 

We have written that last equation in terms of Acj/A/c because this will allow us to label some of 
the velocities as "group" velocities, and others as what I will tentatively call "not group velocities." 
If we have defined ujn = En and /c„(x) = Pn{x), then we can write 

ujn{K,x) = kn-V{x). (5.8) 



If we are considering a semiclassical problem, our weighting coefficients must be well centered on 
some number N with some spread An, such that these satisfy the hierarchy 1 ^ An ^ A^. This 
allows us to define a group velocity for our packet, 



dWr, 



dk-n 



2kN, 



(5.9) 



n=N 



which happens to be the classical velocity of a particle in the potential V with energy E^- If we 
then consider the velocities from Equation ^.7| , writing En = E]\j + and Em = Ej^ + e^, where 
^n,s.m E]\i — V{x) we can simplify our expression, 

{En + Cn) - (En + er. 



^JEn + en- V{x) ±2 yjEN + - V{3 



-■n 



^En - V{x) 1 + 



2[E^-V{x)) (±2(1 + 



2{En-V[x)) 



±l2yjEN - V{x) 



2{En-V{x)){1±2 1) + ( 



±l2^^Ar - V{x) = ±i2kN{x), (-2) 



(5.10) 
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So long as both 6^,6^ ^ E]\r — V{x), as they should be in the semiclassical approximation, then 
half of the velocities contributed by a particular (n, m) will be comparable to the group velocity, Vgr, 
and half will be smaller. It is important to note (in preparation for the next section) that Equation 



5.10| is approximate - it will not happen that half our velocity terms will be degenerate. 

These velocities describe paths (the "traces" in the chapter title) along which the phase is con- 
stant for a given intermode term. That is, given a particular Lnm and a value of x, we can find the 
trajectory along which that phase does not change, and we can do this for all values of x. Although 
what we have done may bear a superficial resemblance to the method of stationary phase, it is very 
different. We could use the method of stationary phase if each intermode term was highly oscillatory 
except along a particular path - we would then assume that each term only made a significant con- 
tribution to the probability density in the immediate vicinity of the path of stationary phase. This 
is invalid here because, typically, the phase term in Equation |5j^ will change at a rate determined 
by the Pn(x), which shouldn't be particularly high. The picture, rather, is of the phase function at 
t = sliding around (and stretching a bit) as time increases. 

5.4 Groups of Velocities and Degeneracy^ 



The result in Equation 5.10| demonstrates that in any problem we will find a range of velocities. 



Those results are, however, approximate. We will show that some intermode terms have exactly 
the same maximum velocities, and that because of that we can treat them as essentially moving 
together. We consider this to be a sort of degeneracy - the sort of degeneracy that produces quantum 
carpets. 



2r 



To this point, my exposition has parlleled that of [Kaplan et al., 200C]. Here, however, I must depart, as their 



treatment of the problem was incorrect. They blithely assert that to discover the traces, we just need to integrate 
from the classical turning points. A cursory survey of the figures in this thesis should be enough to convince you 
otherwise. 
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If we examine the velocities (Equation 5J) at a pointj^ where V{x) = 



±1 



En — E, 



'm 



(5.11) 



nm 




we can factor the numerator 



Vnm — il 



{Ve^. + VEm) {VEi 




Em) 



(5.12) 



and arrive at an important condition, 




(5.13) 



Terms that have the same Vnm at V{x) = have the same classical period - these are our degenerate 
terms. 

We immediately gain some insight into the role of quadratic spectra in producing quantum 
carpets. If the spectrum depends on the the quantum number n squared, then we will have many 
degenerate velocities, while if the spectrum is linear in the quantum number, only traces which 
involve two perfect squares may be degenerate. This suggests that if we begin with a system like 
the simple harmonic oscillator, with a spectrum linear in the quantum number, and set to zero 
all weighting coefficients that are not perfect squares (selecting only states 1,4, 9, 16, etc.), we can 
produce a carpet because we have effectively quadratized^ the spectrum. For an example of this, 
see Figure 5J.. 

Of course, the production of the characteristic canals and ridges of a quantum carpet depends 

not only on degenerate velocities existing, but their producing some sort of peak or valley in the 

^Remember that if we can't find a point where V{x) = 0, we can exploit our ability to change the potential by a 
constant and create one. 

''"Effectively quadratize the spectrum" is my shorthand for choosing weighting coefficients such that only terms 
that have quantum numbers that are perfect squares are contained in the wavefunction. When we do this, we could 
rewrite the spectrum, eigenfunctions, etc., as if they were governed by a new variable, — n, as if they had quadratic 
spectra. Of course, we're not changing the spectrum itself, we're changing the wavefunction. 
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Figure 5.1: Two plots of the simple harmonic oscillator. On the left, the distribution of coefficients 
is Gaussian, with n = 6 and an = 2, while the plot on the right is an even weighting of the perfect 
squares between 1 and 81. 
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probability density. Recalling the form for an intermode term, 

inm{x,t) = ^dnm—^====exp±ii ( [ {±2Pn {x') ±3 Pm {x')) dx' + (En - Em) t) , (5.14) 
^ VPn{x)Pm{x) \J J 

we can see how a sum of these terms might produce a peak or valley at t = 0. First, we ignore 
the imaginary component of the exponential - since our probability density has been reduced to a 
sum of these intermode terms, and the probability density must be real, the imaginary terms will 
ultimately cancel. The sum of intermode terms looks like 



^inra{x,t = ^) = \'S^dnm—j=======CO?>±\\\ {±2Pn{x) ±-iPm,{x')) dx\ , (5.15) 

^ \/Pn\x)Pm.{x) \J J 

or roughly a weighted sum of cosines. Although this hardly constitutes a proolQ that peaks and 
valleys exist, it is not hard to imagine that in many cases we would find them. 

A general calculation of the location of the peaks of an intermode trace would be so general 
as to be useless, and specific examples are likely to be too specific. Instead, we must rest satisfied 
with the notion that very degenerate velocities will typically have enough terms to form some sort 
of peak. As we evolve in time, the peak should roughly follow a classical trajectory, experiencing a 
bit of dispersion. The formula for such a trajectory is 

t{x) = [ V^c/x', (5.16) 

where xq is the location of the center of the peak, and Vnm is some suitable velocity, perhaps the 
middlemost of the degenerate velocities. These are the intermode traces that we are interested in. 
If this is all seems a bit hazy, bear in mind that this idea of a "channel" or "ridge" is something we 
impose on the system, not so easily defined as a local minimum or maximum. We are, in a sense, 
asking a question of the system that it does not want to answer, and this sort of ambiguity is one of 



^Of course, proofs are possible here. I suspect, though, that figuring out what one can prove will be difficult work 
in and of itself. 
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its methods of resistance. Hopefully, future work will uncover a procedure for identifying the most 
significant traces in the carpet. 

5.4.1 A Quick Example 

Once again, we will work on the infinite square well. From the wavefunction. 



^ (x, t) = V2^Cn^ (e+^""« - e-^""M e"*" " *, 



(5.17) 



we can immediately identify the velocities in question as 



Vnm = ±ivr (n ±2 m) 



(5.18) 



We can also find that the wavenumbers k^m will be 



Km = =Fi7r (n =F2 m) 



(5.19) 



so that in any particular velocity bundle there will be a variety of wavenumbers, and we have the 



hope, at least, of interesting interference - for an example, see Figure |5.2| . The most degenerate 



Vnm =P 




0.2 0.4 0.6 ^O.X \ 1 




Vnm=2p 




X 



Figure 5.2: The \vnm\ = tt (left) and \vnm\ = 27r (right) velocity bundles for an even distribution of 
weighting coefficients between n = 1 and n = 10 in the infinite square well. 



velocity will typically be Vnm = i^i", and in time Tr = 2/7r, a trajectory with this velocity will cover 
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a distance of 2 - that is, one full period, so the most prominent traces should have the same period 
as the revival time. Better still, we can ask what velocity we need to have in order to have a period 
equal to the classical period, T^. The condition is 

niT 

Vnm = 27rn, (5.20) 



which should be satisfied by very few (n, m) pairs. An example of how well our separation of the 
wavefunction works is shown in Figures |5.3| - ^^ . 




Figure 5.3: The full carpet, coefficients evenly weighted between n = 1 and n = 10. 
Figure 5.4: The Vnm = intermode terms. Note that they do, indeed, have a period T = Tr. 
Figure 5.5: The Vnm = 27r intermode terms. Note that they have a period T = Tii/2. 



5.5 Degeneracy and ^ci 



Recall that is the "classicized" wavefunction from Section |2J| (see Equation p.21|) , and is what 
governs the spatial distribution and shape of fractional revivals. We have seen that in several cases, 
^ci looks like the original wavepacket moved along a classical path with a relatively small amount 
of dispersion, but we can make this understanding more precise by the application of intermode 
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trace methods. 

Observe from the derivation of Equation ^ that the — Em term comes from the time- 
evolution exponential and the square root terms come from the WKB approximation. If, then, we 
want to find a similar formula for ^ci (see Section 2J and Equation ^.211) , we replace the £"„ — Em 



in the numerator with {2-k/Ti) (n — m). Considering a point where V{x) = 0, we find the following 
velocity degeneracy condition: 

27r n — m 

Vnm = ±17f- ^ , (5.21) 



Unlike Equation |5.13| , we have no hope of factoring this in general. We can, however, see what 
happens in a few obvious cases. If 

En = a^n, (5.22) 



then we quickly find 



2tt n — m , ^ 

±1- 5.23 



a degeneracy condition identical to the one we would have found for the harmonic oscillator in the 
previous section. If, however, we try a spectrum E^ = o?"r? ^ we find something far more interesting. 



27r n — m 

'Jnm = ^ITf; 1 

Iia n ±2 rn 



±i|^, (-2), 

_|_ 2tt n—m / _i_ \ 
1 Tia n+m ' v~'~2J 



(5.24) 



Half of all of the traces are degenerate! Looking at the other half, we want to find two pairs of 
number, (n, m) and {p, q) that will be degenerate. The conditions are 



n-m p-qn + m^O, p + qj^O, 

(5.25) 



n + m p + q 



n ^ p, m ^ p, 
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(p + q) {n — m) = (n + m) {p — q) , 

qn — mp = —qn + mp, 

n p , , 

— = -. 5.26 
m q 

While this is not a particularly difficult equation to satisfy with all of the integers at our disposal, 
it is difficult to satisfy when our weighting coefficients are all within some An ^ n of n » 1 . The 
solutions that will lie closest to (n, m) are {p, q) = 1/2 (n, m) and {p, q) = 2 (n, m). If n and m are 
in the vicinity of n, as they must be in the semiclassical case, that puts p and q near either n/2 or 
2n, both of which would typically be "out of range" of An. For the case of a quadratic spectrum, 
then, we find that half of the velocities are degenerate, and half tend to be non-degenerate. This 
nicely corresponds with our picture of ^ci wavepackets evolving along classical paths (see Figures 



2^ and ^). 

The unfortunate thing about this result for the ^ci of a quadratic spectrum is that it is so 
heavily dependent on the quadraticity of the spectrum. This behavior seems to occur in poten- 
tials with non-quadratic spectra as well, as demonstrated in Figure |5.6| . My short response to 
this is that, in the semiclassical limit, their spectra can appear to be no more than quadratic 
|Nieto and Simmons, 1979|] . This leaves low -n cases to consider, and completing our understanding 



of is a possible subject for future work. 



5.6 Summary 

In this chapter we saw how an analysis of the intermode terms in the probability density of a 
semiclassical wavefunction can give us insight into the degeneracy conditions which must be met in 
order to produce a carpet, and into the carpet itself. When those degeneracy conditions are met, 
we can examine the initial form of certain bundles of velocities and from that find the channels and 
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one with a quadratic spectrum and one with a non 

12 „„x2 



an) 



and on the 

2 



Figure 5.6: Plots of in two potentials 

quadratic spectrum. On the left is the Morse Potential, with En = — {A 

In both 

cases A, B, and a are independent constants. The black lines overlaying the maxima are classical 
paths for those potentials with E = En. For more on these potentials, see Sections A. 2.1 and A.2.6| . 
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ridges of the quantum carpet. For a problem with a spectrum hnear in the quantum number it 
was almost impossible to generate a carpet, though a perverse choice of weighting coefficients could 
effectively quadratize the spectrum. Similar reasoning suggested why from Chapter ^, should 
so often resemble a classically oscillating packet. 



Chapter 6 



Conclusion 

6.1 Results 

The principal result of this thesis may be an understanding of how the different elements in a 
particular problem - potential, eigenfunctions, spectrum, and weighting coefficients - will govern 
revivals, dephasing, and carpet features. The spectrum is the chief determinant of revival dynamics, 
though the initial wavefunction and the potential itself govern the spatial distribution of fractional 
revivals. Though the specifics of dephasing obviously depend on every parameter available, a gross 
dependence on 1/An, where An is the spread in weighting coefficients, seemed to be quite general. 
As for carpet phenomena, we could conclude that the spectrum itself governs the existence of the 
carpet, insofar as it determines whether the velocity degeneracy condition (Equation ^.13|) will be 
met. The weighting coefficients govern (roughly) which channels and ridges are most pronounced, 
meaning that if we were to use the same weighting in two iso-spectral problems, we would expect 
the same number of lines and distribution of line weights. We would not, however, expect the canals 
and ridges to have the same spatial distribution. 
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6.2 Prospects for Future Work 

Although these results are interesting, they are best seen as a platform from which to launch further 
work. 

Are Carpets only a Semiclassical Phenomenon? 

It would appear, from the analysis presented here, that revivals and dephasing are essentially semi- 
classical phenomena, and that carpets are definitely semiclassical phenomena. This is in part because 
our analysis has been directed towards those problems from the beginning. It would be interesting 
to answer conclusively whether a superposition that included a significant weighting of low-energy 
states could produce any of these effects, or demonstrate similar but distinct phenomena. We would 
start the hunt with An « rl 3> 1. 

More than One Quantum Number 

A typical result of increasing the number of quantum numbers in a problem (typically a result of 
increasing the spatial dimensionality) is a proliferation of energy-degenerate states. It is unclear 
whether those new energy degeneracies would contribute to the velocity degeneracy that weaves a 
carpet, or whether their contribution would be unnoticeable. Particularly interesting would be study 
of the two- or three-dimensional harmonic oscillator, which would introduce rather more velocity 
degeneracy than we had before. 

Beyond Quadratic Potentials 

For reasons stated several times, the most important of which is the at-most-quadratic nature of 
the spectrum high in a well, this analysis has been focused on quadratic potentials. It would be 
interesting, though, to see how much of our analysis, and how many of the phenomena, carried over 
into other potentials. This is related to the semiclassical problem - we have many times exploited 
the fact that the limiting spectrum of any well can't rise any faster than as a way of avoiding 
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possibly more interesting (non-n^) behavior near the bottom of the well. 



85 



Open Systems 

We have also only considered bound states in this analysis. We would expect that revival phenomena 
would disappear if probability wasn't being "reflected back" by the walls of the well, but we might 
see revival-like phenomena for a particle that was escaping its well on a time scale longer than its 
revival time. Since they appear relatively early in the time-evolution of the particle, we might also 
see, and perhaps subsequently lose, carpet features. 

Fractals, Fractals, Fractals 

As I stated in the introduction, I still believe that for those wise in the ways of analysis there is 
much to be learned from the self- similarity of quantum carpets. The way is far from clear here 
- perhaps there is a way to go from potential and weighting coefficients to an iterated function 
system? Measure of the fractal dimension of carpets that are not the square well would be a good 
start. 

Other Equations 

I've heard rumor that the Schrodinger equation isn't the only game in town, so far as quantum 
mechanics goes - for example, the Klein-Gordon equation and the non-linear Schrodinger equation. 
What sorts of revival phenomena and carpets do these other equations produce? Do the techniques 
outlined here work there? Related to this, I think, is the application of carpets to BEG. Even if 
they're not related, the existing papers suggest that this is a topic ready for exploration. 

Of course, these techniques may be applied to other wave equations as well. We earlier mentioned 
the correspondence between these problems and waveguides in EM. 
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Finding Particular Intermode Traces 

One prominent problem in Chapter ^ is the lack of a technique to systematically find the locations 
of the peaks and valleys that will be expressed in the carpet. This is not an easy problem, and it 
would be an interesting one to resolve. I suspect that a background in Fourier analysis would be 
helpful in such an attempt. 

Wavepacket Engineering 

There has already been a proposal to use our understanding of revival dynamics for the construction 
of wavepackets with very particular behaviors, | Chen and Yeazell, 1998| . It is not inconceivable that 



we could take this a step further and, using whatever technique we had to find the exact distribution 
of peaks and valleys in a carpet, start from the requirement that a wavepacket have certain canals and 
ridges and come up with a potential and weighting conditions that would produce those features. 
There does not seem to be, as yet, a burning demand for this, but this does sound vaguely like 
something that would be of interest in nanofabrication processes. 



Appendix A 

A Few Exactly Solvable Potentials 



In my numerical work on this project I have rehed heavily on the following exactly solvable poten- 
tials. Knowing the eigenfunctions exactly frees us from the worry that some interesting phenomenon 
is actually an artifact of our numerical integration of the Schrodinger equation. They allow us to 
study a variety of systems - harmomic oscillators, anharmonic oscillators, symmetric potentials, 
non-symmetric potentials, infinite wells, finite wells, double wells, and iso-spectral potentials. These 
have allowed us to test the impact of various properties of the potential on the dynamics - for ex- 
ample, the observation that iso-spectral potentials could have different spatial symmetries ruled out 
the idea that one spectrum (or even one spectrum and one set of weighting coefficients) produced 
one carpet. 

In this Appendix, as in the rest of my thesis, I have set h = 2m = 1. 

A.l Old Friends 



These two potentials are textbook standards - I took them from [[Griffiths, 1995 , 24-44]. I present 
them here to make my notation clear. 
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APPENDIX A. A FEW EXACTLY SOLVABLE POTENTIALS 
A. 1.1 The Infinite Square Well 

The infinite square potential is defined as 



V{x) = < 



0, < X < L, 



oo, otherwise. 



It has a spectrum 



and eigefunctions 



En- ^2 ' 



/2 . (TITX \ 



A. 1.2 The Simple Harmonic Oscillator 

The simple harmonic oscillator potential is 



If we define a dimensionless replacement for x, 



we have a spectrum 



and eigenfunctions 



En= {n+^]u 



V'n(a;) = (^)'^'(2-n!)-i/2^n(0e- 
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A. 2 Cousins from the Old Country 
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These potentials are all taken from a text on supersymmetric quantum mechanics, HCooper et al., 2001 
40-41]. They are all, in supersymmetric-parlance, shape invariant potentials - when subjected to 
the supersymmetry transformation they retain their algebraic form. Many of them make use of the 
following definitions, 

si = s — n + a, 

52 = s — n — a, 

53 = a — n — s, 

54 = — (s + n + a) . 



Each also involves three constants, A, B, and a, along with a few auxiliary functions and constants. 
One warning about these potentials - they are not normalized, though they are normalizable. 
What's more, some of them require a bit of numeric finesse - in Mathematica 4.0, at least, they 
begin to oscillate rapidly in regions where they should be falling off exponentially. 

A. 2.1 The Morse Oscillator 

The Morse oscillator potential is 

V{x) = A^ + B^ exp (-2ax) -2b(^A + ^) exp {-ax) (A.8) 



and has spectrum 



E„ = A"^ - {A- naf 



(A.9) 
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If we define the auxiliary function and constant 



2 B 

y = — e""^, (A.IO) 
a 

s = -, (A.ll) 
a 

we can write the eigenfunctions as 

V'n(x) = y^-"exp (^^ly) L2-2n . ^A.12) 

The Morse oscillator is iso-spectral with the Poschl- Teller potential and the Scarf II (hyperbolic) 
potential. 

A. 2.2 The Eckart Potential 

The Eckart potental is 

V{r) = -\--^-2B coth (ar) + A{A-a) cosech^ar (A. 13) 

and has spectrum 

Er, = A^-{A + naf - f + ^. (A.14) 

[A + na) ^ 

We define the auxiliary function and constants 



y = 


coth (ar 




A 


s = 






a 


A = 


B 






a 




A 


a = 






n + s' 



(A.15) 
(A.16) 
(A.17) 
(A.18) 
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and write the eigenfunctions as 

Mr) = {y- ir^' {y + 1)'^/' Pt'^'Hy)- (a.19) 

A. 2.3 The Poschl- Teller Potential 

The Poschl- Teller potental is 

V{x) = A'^ + [B^ + A^ + Aa) cosech^ (ar) -B{2A + a) coth (ar) cosech (ar) (A. 20) 

and has spectrum 

E^ = A'^ -{A- naf . (A.21) 
We define the auxiliary function and constants 

y = cosh (ar) , (A. 22) 

s = -, (A.23) 
a 

X = (A.24) 
a 

and write the eigenfunctions as 

M^) = {y- 1)^'-'^^' (y + 1)-^'^'^^' p(A-.-i/2,-A -3-1/2) _ (A.25) 

The Poschl- Teller is iso-spectral with the Morse oscillator and the Scarf II (hyperbolic) potential. 

A. 2. 4 The Scarf I (Trigonometric) Potential 

The Scarf I (trigonometric) potential is 

V{x) = -A^ + {A^ + B'^ - Aa) sec^ ax - B {2A - a) tan (ax) sec {ax) (A. 26) 
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and has spectrum 

En = {A + naf + A'^. (A.27) 
Defining the auxihary function and constants 

y = sin ax, (A. 28) 

s = -, A.29 
a 

A = — , (A.30) 
a 

we write the eigenfunctions as 

A. 2.5 The Scarf II (Hyperbolic) Potential 

The Scarf II (hyperbolic) potential is 

V{x) =A^+ {B^ + A^- Aa) sech^aa; + 2B {2A + a) sech (ax) tanh (ax) (A. 32) 

and has spectrum 

E^ = A'^ -{A- naf . (A.33) 
Defining the auxiliary function and constants, 

y = sinhax, (A. 34) 

s = -, A.35 

a 

A = -, (A.36) 
a 
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we write the eigenfunctions as 

^„ = (1 + y^) "'/^ exp (-A tan-^ y) p(iX-s-i/2-iX-s-i/2) _ ^^ 37^ 

The Rosen-Morse II (hyperbohc) potential is iso-spectral with the Poschl- Teller potential and 
the Morse oscillator. 

A. 2. 6 The Rosen-Morse I (Trigonometric) Potential 

The Rosen-Morse I (trigonometric) potential is 

V{x) = A{A-a) cosec^ax + 2B cot ax - A'^ + (A. 38) 



and has spectrum 

— t A J- _ 42 

{A + naf 



q2 ^2 

En = (A + naf - A^ - — To + -f^- (A.39) 



Defining the auxiliary function and constants, 

y = z cot ax, (A. 40) 

s = -, A.41 

a 

A = 4' (A-42) 

a = , (A.43) 

s + n ^ ' 

we can write the eigenfunctions as 
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A. 2. 7 The Rosen-Morse II (Hyperbolic) Potential 

The Rosen-Morse II (hyperbolic) potential is 



V{x) = -A{A- a) cosech^ax + 2B cot ax + A^ + 



and has spectrum 

^2 52 



En = -{A + naf + A^ 



{A + naf ' 
Defining the auxiliary function and constants, 

y = coth ax, 
A 

s = —, 

a 

X - B 

A — —J, 

A 



s + n 



we can write the eigenfunctions as 
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Farey Sequences 



In a typical revival problem, we will only be able to resolve a certain number of revivals, due to 
finite packet width. Because of this the pattern of revivals should follow a Farey sequence. The 
Farey sequence[| is a simple object: F„ is defined as the sequence (in increasing order) of rational 
fractions p/q such that p < q, p and q relatively prime, and q < n. A few examples. 



112 1 

T'3' 2'3'T' 
1 1 2 1 3 2 3 1 

1'4'3' 5' 2' 5' 3'4'T' 
01121323143525341 



(B.l) 
(B.2) 



1'5'4'7'3'8'5'7'2'7'5'8'3'7'4'5'1' ^^"^^ 



If we define Farey addition as 



p n p + n 

-eF- = — , B.4 

q m q + m 



then starting with one Farey sequence we can generate another Farey sequence by Farey adding 
each neighboring term in the given Farey sequence. The resulting fractions will already be reduced, 
and positioned correctly in the sequence. 

In Figure [B.l| we provide a graph of the first 15 Farey sequences, along with a carpet from the 
square well and a quantum beats plot. 



^For more on Farey sequences, see | Harter, 2001 1 or | Weisstein, 1999 1 
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Figure B.l: A plot of the infinite square well (top), with an initial spatial Gaussian wavepacket, 
^ = 1/2, = 0.003, the Farey sequences Fi through F15 (middle), and a two-time-scale quantum 
beat problem with T2/T1 = 200, and a Gaussian distribution of Pk with = 8 (bottom). 
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1.1 A particle in the infinite square well. Lighter areas indicate greater probability den- 



sity. Our initial wavefunction is a Gaussian packet, centered on x = L/2, with) 



ax = 0.003L. For the definition of Tr, see Chapter 2| . 



2.1 II ^11 for a particle in an infinite square well, such that = SOTd. The distribution 



of coefficients is Gaussian, with n = 40 and = 2. For the form of the square well 



that I am using, see Section A.l.l| 20 



2.2 ll^ll for a simple harmonic oscillator. The distribution of coefficients is Gaussian 



with n = 6 and 0"^ = 2. For the version of the simple harmonic oscillator that I am 



usmg, see Section A.lT^ 20 



2.3 Plots of II ^11 for a particle in an infinite square well at t = 0, t = Tr/2, and 



t = Tr/4. Note that the scale is different for the rightmost plot - though conservation 



of probability shrinks our peaks, they have the same shape as the peaks in the other 



two plots. 20 



2.4 Fractional revivals in the infinite square well. Our initial wavefunction is a Gaussian 



packet, centered on a; = -L/2, with ax = 0.003L. The symmetry of our wavepacket 



allows it to have a full revival (modulo an overall phase factor) at t = Tr/8. At times 



below that, though, we should be able to count the fractional revivalsj 22 
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2.5 V{ar) (above) and ^ci ^ (below) for a particle in the Poshcl- Teller potential (see| 
Appendix Aj ). The dotted line in the above plot indicates En, in the lower plot they 



indicate four different classical paths of energy En. The Poschl- Teller potential has 



a purely quadratic spectrum. For more details on it, see Section |A.2.g 



2.6 V(ax) (above) and H^efll (below) for a particle in the Morse potential (see Appendix 
A| ). The dotted line in the above plot indicates En, in the lower plot they indicate four 



different classical paths of energy En. The Morse potential has a purely quadratic 



spectrum. For more details on it, see Section |A.2.1 



2.7 for roughly Gaussian spatial distribution, centered at x = a/3 



3.1 A Gaussian distribution of the Pk, with T2/T1 = 200 and (j„ = 8. In the plot of the 



early evolution, by Equation 3.2E dephasing should occur at t ^ 3.125ri 



3.2 A Gaussian distribution of the Pj., with T2/T1 = 200 and N = 8. In the plot of the 
early evolution, by Equation 3^ dephasing should occur at t 3.125Ti 

3.3 Re-{Erf(0x)[, X g 1-10, 10| 



3.4 Im -{Erf ((/>x)}, X £ 1-10, 10|. 



4.1 These density plots show a solution of the infinite square well in the region r £ 



[0,4], g £ [0,1] for An = 5, n = 30. On the left are the five relevant terms from 



Equation 4.27| , on the right are the middle 10 terms of the sum in Equation |4.2 



4.2 All of these show the region r g [0.5, 1.5], ^ £ [0, 1] for An = 5, n = 30, and the time 



it took Mathematica to plot them. On the left is the I = —1 term from Equation 



4.27, the middle figure contains the middle 20 terms of the sum in Equation 4.2, and 



at right is a plot of the middle 10 terms from that same Equation. 



/ \ 

4.3 A plot of exp —7^ (dashed) and sech jaT^) (solid), —5a < y < 5a. The constant 



a = arcsech (-^) ^ 2.69 is chosen so that the two functions are equal at y = ±2a. 
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4.4 These density plots show a solution of the infinite square well in the region r £ 



0, 4J, g £ [0, IJ for An = 5, n = 30, and the white line indicates the lower of the 



A = 1 solutions of Equation 4.45|. On the left are the five relevant terms from 



Equation |4.27| , on the right are the middle 10 terms of the sum in Equation |4.2| . Note 
that the ridges in the left plot only appear above the line, and in the right plot that 



this holds for ^ < 1/2 - for ^ > 1/2, the interference between terms in Equation 4.27 
becomes non-negligible 66 



5.1 Two plots of the simple harmonic oscillator. On the left, the distribution of coefficient; 



is Gaussian, with re = 6 and (7^ = 2, while the plot on the right is an even weighting 



of the perfect squares between 1 and 81. 75 



5.2 The v„, 



TT (left) and 



27r (right) velocity bundles for an even distribution 



of weighting coefficients between re = 1 and n = 10 in the infinite square well.| .... 77 



5.3 The full carpet, coefficients evenly weighted between re = 1 and re = 10. 78 



5.4 The Vnm = intermode terms. Note that they do, indeed, have a period T = Tjj 



5.6 Plots of \&ef in two potentials, one with a quadratic spectrum and one with a non- 



7J 



5.5 The Vnm = ^tt intermode terms. Note that they have a period T = Tr/2. 78 



quadratic spectrum. On the left is the Morse Potential, with En = A — (A — an) 



and on the right the Rosen-Morse I potential, with £*„ 



-A-^ + (A + an) 



[Bl {A + an)Y + {B/AY. In both cases A, B, ando; are independent constants 



The black lines overlaying the maxima are classical paths for those potentials with 

E = En. For more on these potentials, see Sections A.2l| and |A.2.6| 81 



B.l A plot of the infinite square well (top), with an initial spatial Gaussian wavepacket 



g = 1/2, = 0-003, the Farey sequences Fi through F15 (middle), and a two-time 
scale quantum beat problem with T2/T1 = 200, and a Gaussian distribution of 1% 
with (7fc = 8 (bottom) 



96 



Bibliography 



[Abramowitz and Stegun, 1965] Abramowitz, M. and Stegun, I. A. (1965). Handbook of mathemat- 
ical functions. Dover. 

[Aronstein, 2000] Aronstein, D. L. (2000). Analytical investigation of revival phenomena in the 
finite square-well potential. Physical Review A, 62:022102. 

[Aronstein and Stroud, 1997] Aronstein, D. L. and Stroud, C. R. J. (1997). Fractional wave-function 
revivals in the infinite square well. Physical Review A, 55(6):4526-4537. 

[Averbukh and Perehnan, 1989] Averbukh, I. S. and Perelman, N. (1989). Fractional revivals: Uni- 
versality in the long-term evolution of quantum wave packets beyond the correspondence principle 
dynamics. Physics Letters A, 139(9) :449-453. 

[Averbukh and Perelman, 1991] Averbukh, I. S. and Perelman, N. (1991). The dynamics of wave 
packets of highly-excited states of atoms and molecules. Sov. Phys. Usp., 34(7):572-591. 

[Berry, 1996] Berry, M. V. (1996). Quantum fractals in boxes. Journal of Physics A: Mathematical 
and General, 29:6617-6629. 

[Berry and Bodenschatz, 1999] Berry, M. V. and Bodenschatz, E. (1999). Caustics, multiply recon- 
structed by talbot interference. Journal of Modern Optics, 46(2):349-365. 

[Berry and Klein, 1996] Berry, M. V. and Klein, S. (1996). Integer, fractional, and fractal talbot 
effects. Journal of Modern Optics, 43(10):2139-2164. 

100 



BIBLIOGRAPHY 101 

[Bluhm et al., 1996] Bluhm, R., Kostelecky, V. A., and Porter, J. A. (1996). The evolution and 
revival structure of quantum wave packets. American Journal of Physics, 64:944. 

[Chen and Yeazell, 1998] Chen, X. and Yeazell, J. A. (1998). Analytical wave-packet design scheme: 
Control of dynamics and creation of exotic wave packets. Physical Review A, 57(4):R2274-R2277. 

[Choi et al., 2001] Choi, S., Burnett, K., Friesch, O. M., Kneer, B., and Schleich, W. (2001). Spa- 
tiotemporal interferometry for trapped atomic bose-einstein condensates. Physical Review A, 
63:065601. 

[Cooper et al., 2001] Cooper, F., Khare, A., and Sukhatme, U. (2001). Supersymmetry in Quantum 
Mechanics. World Scientific Publishing Co. 

[Courant and Hilbert, 1953] Courant, R. and Hilbert, D. (1953). Methods of mathematical physics. 
Interscience Publishers, New York,, 1st english edition. 

[Dubra and Ferrari, 1999] Dubra, A. and Ferrari, J. A. (1999). Diffracted field by an arbitrary 
aperture. American Journal of Physics, 67(l):87-92. 

[Friesch et al., 2000] Friesch, O., Marzoli, I., and Schleich, W. P. (2000). Quantum carpets woven 
by wigner functions. New Journal of Physics, 2(4):4.1-4.11. 

[Griffiths, 1995] Griffiths, D. J. (1995). Introduction to quantum mechanics. Prentice Hall, Engle- 
wood Cliffs, N.J. 

[Grossmann et al., 1997] Grossmann, F., Rost, J.-M., and Schleich, W. P. (1997). Spacetime struc- 
tures in simple quantum systems. Journal of Physics A: Mathematics and General, 30:L277-L283. 

[Hall et al., 2001] Hall, M. J., Reineker, M. S., and Schleich, W. P. (2001). Unravelhng quantum 
carpets: A travelling wave approach. arXiv, (quant-ph/9906107 v2). 

[Barter, 2001] Harter, W. G. (2001). Quantum-fractal revival structure in cn quadratic spectra: 
Base-n quantum computer registers. Physical Review A, 64:012312-1. 



BIBLIOGRAPHY 102 

[Jie et al., 1998] Jie, Q.-L., Wang, S.-J., and Wei, L.-F. (1998). Partial revivals of wave packets: 
An action-angle phase-space description. Physical Review A, 57(5):3262-3267. 

[Kaplan et al., 2000] Kaplan, A., Marzoli, I., Lamb, w. J., and Schleich, W. (2000). Multimode 
interference: Highly regular pattern formation in quantum wave-packet evolution. Physical Review 
A, 61:032101. 

[Kaplan et al., 1998] Kaplan, A., Stifter, P., van Leeuwen, W., Lamb, W. J., and Schleich, W. P. 
(1998). Intermode traces-fundamental interference phenomenon in quantum and wave physics. 
Physica Scripta, T76:93-97. 

[Knospe and Schmidt, 1996] Knospe, O. and Schmidt, R. (1996). Revivals of wave packets: General 
theory and application of rydberg clusters. Physical Review A, 54(2):1154-1160. 

[Leichtle et al., 1996a] Leichtle, A., Averbukh, L S., and Schleich, W. (1996a). Generic structure of 
multilevel quantum beats. Physical Review Letters, 77(19) :3999-4002. 

[Leichtle et al., 1996b] Leichtle, C., Averbukh, L S., and Schleich, W. (1996b). Multilevel quantum 
beats: An analytical approach. Physical Review A, 54(6):5299-5312. 

[Lock and Andrews, 1992] Lock, J. A. and Andrews, J. H. (1992). Optical caustics in natural 
phenomena. American Journal of Physics, 60(5):397-407. 

[Loinaz and Newman, 1999] Loinaz, W. and Newman, T. (1999). Quantum revivals and carpets in 
some exactly solvable systems. Journal of Physics A: Mathematics and General, 32:8889-8895. 

[Marzoli et al., 1998a] Marzoli, L, Bialynicki-Biruli, L, Friesch, O., Kaplan, A., and Schleich, W. 
(1998a). The particle in the box: Intermode traces in the propagator. arXiv, (quant-ph/9804015). 

[Marzoli et al., 1998b] Marzoli, L, Saif, F., Bialynicki-Birula, L, Friesch, O., Kaplan, A., and Schle- 
ich, W. (1998b). Quantum carpets made simple. Acta Physica Slovaca, 48(3):323-333. Wigner 
Functions. 



BIBLIOGRAPHY 103 

[Nieto and Simmons, 1979] Nieto, M. M. and Simmons, L. J. (1979). Limiting spectra from confin- 
ing potentials. American Journal of Physics, 47(7):634-635. 

[Provost and Baranger, 1993] Provost, D. and Baranger, M. (1993). A semiclassical calculation of 
scars for a smooth potential. arXiv, (cliao-dyn/9305003). 

[Razi Naqvi et al., 2001] Razi Naqvi, K., Waldenstrom, S., and Haji Hassan, T. (2001). Fractional 
revival of wave packets in an infinite square well: a fourier perspective. European Journal of 
Physics, 22:395-402. 

[Rozmej and Arvieu, 1998] Rozmej, P. and Arvieu, R. (1998). Clones and other interference effects 
in the evolution of angular-momentum coherent states. Physical Review A, 58(6):4314-4329. 

[Ruostekoski et al., 2001] Ruostekoski, J., Kneer, B., Schleich, W. P., and Rempe, G. (2001). In- 
teference of a bose-einstein condensate in a hard-wall trap: Prom the nonlinear talbot effect to 
the formation of vorticity. Physical Review A, 63:043613. 

[Saif, 2000] Saif, F. (2000). Quantum recurrences: Probe to study quantum chaos. Physical Review 
E, 62(5):6308-6311. 

[Silverman, 1995] Silverman, M. P. (1995). More than one mystery : explorations in quantum 
interference. Springer- Verlag, New York. 

[Weisstein, 1999] Weisstein, E. W. (1999). Farey sequences, 

http:/ /mathworld. wolfram.com/FareySequences. html. 

[Wojcik et al., 2000] Wojcik, D., Bialynicki-Birula, I., and Zyczkowski, K. (2000). Time evolution 
of quantum fractals. Physical Review Letters, 85(24) :5022-5025. 

[Wojcik and Zyczkowski, 2001] Wojcik, D. and Zyczkowski, K. (2001). Fractality of certain quan- 
tum states. arXiv, (math-ph/0107030). 

[Wolfram, 2002] Wolfram (2002). The error function. http://functions.wolfram.com/GammaBetaErf/Erf/. 



BIBLIOGRAPHY 104 

[Wright et al., 1997] Wright, E. M., Wong, T., Collett, M. J., Tan, S. M., and Walls, D. F. (1997). 
Collapses and revivals in the interference between two bose-einstein condensates formed in small 
atomic samples. Physical Review A, 56(l):591-602. 



